Dirac fermions in strong electric field and quantum transport in graphene 
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Our previous results on the nonperturbative calculations of the mean current and of the energy- 
momentum tensor in QED with the T-constant electric field are generalized to arbitrary dimensions. 
The renormalized mean values are found; the vacuum polarization and particle creation contributions 
to these mean values are isolated in the large T-limit, the vacuum polarization contributions being 
related to the one-loop effective Euler-Heisenberg Lagrangian. Peculiarities in odd dimensions are 
considered in detail. We adapt general results obtained in 2 + f dimensions to the conditions 
which are realized in the Dirac model for graphene. We study the quantum electronic and energy 
transport in the graphene at low carrier density and low temperatures when quantum interference 
effects are important. Our description of the quantum transport in the graphene is based on the 
so-called generalized Furry picture in QED where the strong external field is taken into account 
nonperturbatively; this approach is not restricted to a semiclassical approximation for carriers and 
does not use any statistical assumptions inherent in the Boltzmann transport theory. In addition, 
we consider the evolution of the mean electromagnetic field in the graphene, taking into account the 
backreaction of the matter field to the applied external field. We find solutions of the corresponding 
Dirac-Maxwell set of equations and with their help we calculate the effective mean electromagnetic 
field and effective mean values of the current and the energy-momentum tensor. The nonlinear 
and linear I — V characteristics experimentally observed in both low and high mobility graphene 
samples is quite well explained in the framework of the proposed approach, their peculiarities being 
essentially due to the carrier creation from the vacuum by the applied electric field. 
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I. INTRODUCTION 
A. General 

It is well-known that QFT with an external back- 
ground is an adequate model for studying quantum pro- 
cesses in the cases when a part of the quantized field is 
strong enough to be treated as a given and a classical one. 
Numerous problems in QED and QCD with superstrong 
electromagnetic fields, which must be treated nonpertur- 
batively, are at present investigated in this framework, 
with applications to astrophysics and condensed matter 
physics (e.g. graphene physics) (Trill- In these models, 
one often needs to know how the external background af- 
fects admitted states of the quantized fields and to what 
extent the external background idealization is consistent, 
i.e., whether the backreaction on the external background 
is sufficiently small. Here we consider the quantized Dirac 
field interacting with an external electromagnetic back- 
ground. The external electromagnetic field changes the 
states of the Dirac field and induces transitions between 
them. These changes in general occur as vacuum polar- 
ization effects combined with those due to particle cre- 
ation. A uniform constant magnetic field, for instance, 
can produce only vacuum polarization effects, whereas 
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particle creation is due to electric-like external fields. 



The first calculations of vacuum polarization effects 
have been carried out by Heisenberg and Euler in the case 
of QED with constant, uniform, parallel electric E and 
magnetic B fields, with E <C E c = m 2 /e ~ 10 16 V/ cm 
{E c is Schwinger's critical field) and arbitrary B, ignor- 
ing particle creation It turns out that the vacuum 
polarization effects (which are local in time in the uni- 
form background) can be described by some effective 
Lagrangians, in particular by the Heisenberg-Euler La- 
grangian, see @. On the other hand, the particle cre- 
ation effects are global ones, as they depend on the his- 
tory of the external field action and, therefore, cannot 
be described by any local effective Lagrangian (this does 
not contradict the fact that the imaginary part of the 
Heisenberg-Euler effective action is related to the total 
number of particle created, since this action is a global 
quantity). A non-perturbative approach (with respect 
to the external background) to the particle creation in 
the so-called f-electric steps (electric-like fields that are 
switched on and off at initial and final time instants, re- 
spectively) was elaborated in the framework of relativis- 
tic quantum mechanics in Q, and in the framework of 
QFT in [l|; a more complete list ofpublications on the 
particle creation can be found in [3(. A direct way to 
study the influence of an external electromagnetic back- 
ground on the Dirac field is calculating nonperturbatively 
the mean energy-momentum tensor. In case of a quasi- 
constant electric field (the so-called T-constant electric 
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i.e. a uniform electric field, constant within a time 
interval T, and zero outside it) such a problem was solved 
in 0, [l(|. The obtained result allows one to isolate, un- 
der certain conditions, vacuum polarization effects from 
the particle creation effects. 

Until recently, problems related to particle creation 
from the vacuum were of a purely theoretical interest. 
This is related to the fact that, due to the presence of 
large gaps between the upper and lower branches in the 
spectrum of Dirac particles, particle creation effects can 
be observed only in huge external electric fields of the 
magnitude of E c . However, recent technological advances 
in laser science suggest that lasers such as those planned 
for the Extreme Light Infrastructure project (ELI) may 
be able to reach the non-perturbative re gim e of pair pro- 
duction in the near future (see review Moreover, 
the situation has changed completely in recent years re- 
garding applications to condensed matter physics: par- 
ticle creation became an observable effect in graphene 
physics, an area which is currently under intense de- 
velopment [TH, [l|. Briefly, this is explained by two 
facts: first, the low-energy electronic excitations in the 
graphene monolayer in the presence of an external elec- 
tromagnetic field can be described by the Dirac model 
pjj], namely, by a 2 + I quantized Dirac field in such a 
background (that is, dispersion surfaces are the so-called 
Dirac cones); and, second, the gap between the upper 
and lower branches in the corresponding Dirac particle 
spectra is very small, so that the particle creation effect 
turns out to be dominant (under certain conditions) as a 
response to the external electric-like field action on the 
graphene. In particular, this effect is crucial for under- 
standing the conductivity of grap hene, specially in the 
so-called non- linear regime 0, H3-[3- The first experi- 
mental observation of non-linear current- voltage charac- 
teristics (I—V) of graphene devices and its interpretation 
in terms of pair-creation has been recently reported [l9| . 

It was also shown recently in [20] that linear disper- 
sion surfaces with the Dirac points, similar to those in 
the graphene monolayers, are generic in the spectrum of 
non-relativistic 2 + 1-dimensional Hamiltonians with po- 
tentials that have the symmetry of a honeycomb struc- 
ture. Thus, one may expect the discovery of new ma- 
terials in which the conductivity is described by a Dirac 
model of (2 + 1) massless fermions. Moreover, a mass gap 
in the graphene band structure can be generated by sev- 
eral methods. One example is given by graphene nanorib- 
bons. In this case one has a quasi-one-dimensional spatial 
geometry that confines the graphene electrons to a strip 
of large length and small width. The confinement gap de- 
pends on the strip width and on imposed boundary con- 
ditions, see [HI for a review. The spectrum of graphene 
can also be gapped by explicit or spontaneous sublattice 
symmetry breaking, see |21| and references therein. It is 
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an important fundamental and practical problem under 
current research. 

It should also be noted that there is a deep connec- 
tion between graphene, which has two Dirac cones, and 
topological insulators, which are characterized by a sin- 
gle Dirac cone on each surface; see [IH [HJ for a review 
and [23l | for the first experimental realization of a single- 
valley Dirac system in zero-gap HgTe quantum wells. In 
these Dirac systems, the mass gap can be generated by an 
appropriate selection of material parameters that depend 
on the quantum well geometry. In our consideration, it 
is assumed that the two cones of graphene are decoupled 
and the system behaves like two copies of a single Dirac 
cone. Thus, the results obtained for graphene should 
also be relevant for the single Dirac cone on the surface 
of a topological insulator. Then we expect the interface 
transport properties of topological insulator to be similar 
to those described for graphene. 

In the present article, we adapt our general results on 
the particle creation to the cases which are realized in 
the Dirac model for graphene. First, our previous results 
in calculating the vacuum mean current and the energy- 
momentum tensor of the Dirac field in the T-constant 
electric field 0, [l(| is extended to arbitrary dimensions, 
especially to (2+l)-d. Then we consider some of immedi- 
ate applications of these results to a number of problems 
related to the dc conductivity of graphene in the so-called 
superlinear (non-perturbative) regime. 

We show that the application of non-perturbative 
methods of QED with strong field and unstable vac- 
uum to studying the Dirac model in graphene allows 
one to more adequately describe quantum transport in 
graphene. The matter is that usually the electronic trans- 
port in graphene, even near the Dirac point, is described 
within the framework of traditional methods used in con- 
densed matter physics. For example, these are methods 
based on WKB description of the carriers [TH HH , or nu- 
merical methods that attend one-particle quantum me- 
chanical description of the carriers [I?) ■ Other numerical 
methods exploit a Green's function formalism of non- 
equilibrium statistical mechanics adopted for describing 
states not very far from the equilibrium 19]. In some 
cases the Boltzmann transport theory, based on WKB 
approximation, describes the system evolution quite well, 
and sometimes it allows one to find analytic solutions. 
However, these cases are restricted to the ones when the 
large carrier densities obey the assumption that charac- 
teristic length scales of the system be larger than typi- 
cal wave lengths of carriers. Moreover, the Boltzmann 
theory is not adapted to describe motion of massless 
carries created from the vacuum. The numerical meth- 
ods in their present form are not adjusted for study- 
ing the time-evolution of the system that would take 
the electrodynamic backreaction into account, whereas 
the latter may be of essential importance in the present 
highly non-equilibrium environment caused by the pair 
creation. Moreover, the electric field creates carriers in 
pure states, whose distribution differs significantly from 
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the equilibrium distribution. This is a nonstandard situ- 
ation for usual transport problems in the condensed mat- 
ter physics. Thus, a proper description of the quantum 
transport in graphene close to the Dirac point both in 
the ballistic case and in the presence of a disorder is still 
an open problem. 

Our description of the quantum transport in the 
graphene is based on strong-field QED, it is not restricted 
by a semiclassical approximation of carriers and it does 
not use any statistical suppositions typical of the Boltz- 
mann transport theory. Our approach is based on exact 
solutions of the Dirac equation, where the strong external 
field is taken into account non-perturbatively. A strong 
field approximation used for analytical calculations is re- 
lated to the consistency of the Dirac model with a given 
external field and we show that it is well-working under 
certain conditions. In fact, we study the electronic and 
energy transport in graphene at low carrier density and 
low temperatures when quantum interference effects are 
important. 

We consider the evolution of the mean electromagnetic 
field in the graphene, taking into account the backreac- 
tion of the quantum matter field on the applied exter- 
nal field. In doing this, we use consistency restrictions 
that describe the regime when the backreaction can be 
neglected. We derive restrictions from above on the al- 
lowed strength of the external electric field and on its du- 
ration, admitted at its given strength. In making some 
experimental conclusions, we try to compare these re- 
strictions with typical experimental scales. We present 
a generalization of the QED model with the T-constant 
external field to take into account the backreaction of 
the mean current to the applied electric field. We find a 
self-consistent solution of the Dirac-Maxwell set of equa- 
tions for this generalized model and calculate the effec- 
tive mean field and effective mean values of the current 
and the energy-momentum tensor. We show that the 
non-linear and linear I — V experimentally observed in 
low and high mobility samples, respectively, can be ex- 
plained in the framework of the presented consideration, 
and that such a behavior is a consequence of the fact that 
the conductivity in the graphene is essentially due to the 
pair creation from the vacuum by the applied electric 
field. 

The article is organized as follows: In Subsection ll Bl we 
give introductory overview of the basics of the QED with 
unstable vacuum and introduce the necessary general no- 
tation. In Section |TT1 we introduce the Dirac model in a 
T-constant background, we fix notation and collect pre- 
vious results (mainly of [ToL l24j) required for our work. 
Exact solutions of the Dirac equation for this background 
are presented, and we describe how they are related to 
the particle creation and to the mean values of field ob- 
servables, following the general theory presented in Q. In 
Section IIIII we consider the vacuum mean values of the 
current density and the energy-momentum tensor in a 
T-constant background for the case of arbitrary dimen- 
sions. Namely, in Subsection 1111 Al we express the vac- 



uum mean values of the current density and the energy- 
momentum tensor in terms of appropriate Green func- 
tions constructed by the help of the exact solutions of the 
Dirac equation in the external field. In Subsection IIIIBI 
peculiarities in odd dimensions are considered. In Sec- 
tion [nTC] vacuum polarization and particle creation con- 
tributions to the mean values of the energy-momentum 
tensor are isolated in the large duration approximation, 
and the vacuum polarization contributions are related to 
the one-loop effective Euler-Heisenberg Lagrangian. In 
Subsection 1111 Dl we find renormalized mean values of the 
energy-momentum tensor using the zeta-function regu- 
larization and compute the vacuum polarization contri- 
butions to these mean values in arbitrary dimensions. In 
Subsection IIII El we compute the particle creation con- 
tributions to the mean values of the current and of the 
energy-momentum tensor in arbitrary dimensions to gen- 
eralize the four dimension results of Pair-creation 
contributions are finite, due to the finite duration of the 
field. In Section HVl our results are applied to graphene 
physics. Namely, in Subsection IIV Al the Dirac model of 
graphene is briefly described, and the mean values com- 
puted in the previous section are used to determine the 
mean current and energy-momentum tensor in the mate- 
rial. In Subsection llV Bl these results allow us to estimate 
scales for backreaction effects, using the appropriate so- 
lution of the Maxwell equations. We show that the super- 
linearity of I— V observed experimentally in low mobility 
samples is essentially due to the pair creation from the 
vacuum by the applied electric field in the regime where 
the backreaction is negligible. In Subsection IIV CI we 
study a generalization of the model with the T-constant 
external field to taking into account the backreaction of 
the mean current to the applied electric field. We find a 
self-consistent solution of the Dirac-Maxwell set of equa- 
tions (with unstable vacuum in fcrmion sector) and cal- 
culate the effective mean field and effective mean values 
of the current and of energy-momentum tensor. We show 
that linear I — V experimentally observed in high mobil- 
ity samples is due to the pair creation from the vacuum 
by the applied electric field in the regime where backre- 
action is important. In Section [V] (Summary) we briefly 
list the main new results obtained in the article and add 
some relevant comments. 



B. Dirac model with external background 

In order to study the Dirac field in the external back- 
ground, we use a formulation of QED with unstable vac- 
uum (the generalized Furry representation) developed in 
0, [25| , where a strong external field is treated nonpertur- 
batively. In this work we consider the case of the so-called 
T-constant electric field, i.e., we assume that for t < ti n 
and for t > t out , the T-constant electric field is absent, 
therefore initial and final vacua are vacuum state of free 
in- and out- particles respectively. During the interval 
tout —tin = T, the Dirac field interacts with a constant 
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uniform electric field. The initial and final vacua are dif- 
ferent due to the difference of the initial and final value 
of the external electromagnetic field potentials. In the 
Heisenberg picture, there exists a set of creation and an- 
nihilation operators a n (in), a n (in) of in-particles (elec- 
trons), and operators b n (in), b n {in) of m-antiparticlcs 
(positrons), the corresponding in- vacuum being |0, in) at 
the same time there exists a set of creation and annihila- 
tion operators ajj(out), a n (out) of out-electrons and oper- 
ators b n (out), b n (put) of out-positrons, the corresponding 
out- vacuum being |0, out), 

a n (in)\0,in) = b n (in)\0,in) — 0, Vn, 
a n (out)\0, out) = b n (out)\0,out) = 0, Vn. 

In both cases, by n we denote complete sets of quantum 
numbers describing in- and out- particles. The in- and 
out-operators obey the canonical anticommutation rela- 
tions: 

[a n (in),a n ,(in)] + = [a n (out),a n ,(out)] + = 8 n . n , , 
[b n (in),b' n ,(in)} + = [b n (out),b n ,(out)}+ = 8 n ^ n > , 

the remaining anticommutators being zero. 

The m-operators are associated with a complete or- 
thonormal set of solutions {t;tpn{x)} (C = + f° r electrons 
and £ = — for positrons) of the Dirac equation with In- 
constant electric field. Their asymptotics as t < ti n can 
be classified as free particles and antiparticles. The out- 
operators are associated with a complete orthonormal 
out-set of solutions {''ipn (%)} of the Dirac equation with 
T-constant electric field. Their asymptotics as t > t out 
can be classified as free particles and antiparticles. The 
in- and out- operators are defined by the two representa- 
tions of the quantum Dirac field ^(x) in the Heisenberg 
representation 

^f(x) = 2J [a n (in) +ip n (x) + b^in) -i/j n (x)] 

n 

= ^ [a n (out) + ijj n (x) + b'l(out)) ~tp n (x)] . (1) 

n 

In- and out-solutions with given quantum numbers n are 
related by a linear transformation of the form: 

^ n (x)=g( + \<)+iP n (x)+g(^)-i; n (x) , (2) 

where the g's are some complex coefficients. Then a 
linear canonical transformation (Bogolyubov transforma- 
tion) between in- and out- operators which follows from 
Eq. ([TJ is defined by these coefficients. 

The vacuum mean current vector, energy and momen- 
tum vector of the quantum Dirac field Vt'(x) at a time 
instant t are defined as integrals over the spatial volume. 
Due to translational invariance of the external field under 
consideration, all these mean values are proportional to 
the space volume. Therefore, it is enough to calculate the 
vacuum mean values of the current density vector (y (t)) 



and of the energy- momentum tensor (EMT) (T^^t)), 

(y(t)) = (Q,in\f%in), 

(T M „(i)> = (0,in\T^\0,in) . (3) 

Here we stress the time-dependence of mean values (|3|), 
which does exist due to the time-dependence of the ex- 
ternal field. The operators of the current density and the 
EMT are described in terms of the quantum Dirac field 
as follow^ 

f = I [*(*), y *(!)] , T M „ = \(T™ + T-") , 

T™ n = \ {[* (x),^P^(x)} + [P;#(z) )7 ^0r)]} , 
P„ = ify - qA^(x), *(s) = * t (a;)7 , 

where A^(x) are electromagnetic potentials of the exter- 
nal field, q is the particle charge (for an electron q = — e), 
and ^(x) is the Heisenberg operator of the Dirac field 
that obeys the Dirac equation with the external back- 
ground. 

Note that the mean values ^ depend on the definition 
of the initial vacuum, \0,in) and on the evolution of the 
electric field from the time ti n of switching it on up to 
the current time instant t, but they have nothing to do 
with the further history of the system. The renormal- 
ized vacuum mean values (y (t)) and (T^„(t)), U n < t < 
t out are sources in equations of motion for mean electro- 
magnetic and metric fields, respectively. In particular, 
complete description of the backreaction is related to the 
calculation of these mean values for any t. 

II. DIRAC FIELD IN T-CONSTANT ELECTRIC 
BACKGROUND 

In what follows we are going to deal with the Dirac 
equation in (d = D + I)-dimensional Minkowski space 
with an external electromagnetic field A fi (x), 

(7"P M - m)ip(x) = . (4) 

Here ip(x) is a 2 I d / 2 1 -component spinor (the brackets 
stand for 'integer part of), m is the mass, q is the charge, 
the Greek index assumes values u = 0,1, ... ,D, and the 
gamma matrices satisfy the standard anti-commutation 
relations: 

[y , 7*1+ = 2rT, ti mv = diag(l, -1, . . . , -1). 

In what follows, we consider the T-constant field de- 
scribed by a vector potential with only one nonzero com- 
ponent Ai(t) (A^t) = 0, u ^ I), 

{tin £ £ I ( oo, £^ n ), ti n — T /2 , 
t, t€lnt = [t in , t ou t], 
touU t G II = {t ou t, °°) i tout = T/2 . 



2 We use the relativistic units h = c = 1, in which the fine structure 
constant is a = e 2 /ch = e 2 . 
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such that the electric field E (t) has also only one nonzero 
component, which is nonzero for t £ Int, i.e., 

E 1 (t) = E, t e Int; E l (t) = 0, t € IUII. 

If one represents the spinor ip(x) in the form 

iP(x) = (7^ + m)0(x) , (5) 

where <f>(x) is a new spinor, then <f>{x) obeys the following 
equation: 



where £ = ±, 



iA/2, and 



P 2 - to 2 - 



(6) 



where 



= c^A, - = -[y, y] , 

|a^P^ - i9 P(i) 7 V . 

Solutions of the Dirac equation in a T-constant field were 
studied in detail in (24|. Below, we use these results. 

First we choose a set of constant orthonormalized 
spinors v a , r , with s = ±1, and r = (n, f2, . . . , r[d/2]-i), 
Ti = ±1, such that 7 7 1 w S!r = su Sir . The indices de- 
scribe the spin polarization, which is not coupled to the 
electric field, and together with the additional index s 
provide a suitable parametrization of the solutions. Then 
we represent the spinors <p(x) as follows: 



(7) 



Thus, the time-evolution is described by the ordinary dif- 
ferential equation of second order 

+ ]pi ~ qAi(t)} 2 + \qE\ A + isqE(t)} <p p , s (t) = , 

(8) 

where A = (pj_ + m 2 )/\qE\ and pj^ is the transversal 
momentum, p^ = (0,p 2 , . . . ,p D ). 

At early (t < t in -region I), and late (t > t out -region 
II) times, Eq. (|5]) has plane wave solutions £</> PlS ,r and 
^0p,s,r respectively, with £ = ±1, which satisfy simple 
dispersion relations 

I: c^,-(*)~e-« w *»*, II: € <p v ,.(t) ~ e^"""* , 

u m / ou t = \J{pi - qEt in/out ) 2 + p 2 ± + m 2 . (9) 

As was demonstrated in [24j . the in-set {^J p . r } and out- 
set {^ip p , r } °f solutions of the Dirac equation in the In- 
constant electric field can be taken in the form 

±4> P ,r(x) = (7 • P + m)±0p,±, r (a:) , 

± V Pl r(^) = (7 • P + f») ± ^ PjT , r (a:) . (10) 

For t £ Int, the general solution of Eq. © is completely 
determined by an appropriate pair of the following linear 
independent Weber parabolic cylinder functions: 

AMl+«)/2 (CI 1 - 00 > I>_v-(l-.)/2 (C(i + 00 , 



|g£|f-pisgn (gg) 



According to relation ([2]), an owt-solution corresponding 
to a plane wave in interval II is thus described by a su- 
perposition of the Weber functions in interval Int, and 
extends into a superposition of particle and antiparticle 
m-solutions in I. Then one can explicitly find the coef- 
ficients g. In terms of such coefficients, the differential 
mean number H p r of particles created from vacuum with 
given p and r at a time instant t > t ou t is 

ttp.r = (Q,in\a P r (out)a Pir (out)\0,in) = N p = |.g(-| + )| 2 , 

where the Bogolyubov transformation between in- and 
out- creation and annihilation operators is used. The 
number of particles created is equal to the numbers of 
antiparticles created. Then K Pjr can be treated as the 
number of pairs created. The result in (jlip is indepen- 
dent of the spin polarization. That is why we use in what 
follows the notation H p for the quantity (fTTj) . Note that 
there is no particle production after the time instant t ou t- 
Thus, H p depends only on the interval T. The explicit 
expression for K p was studied in detail in [2J|. Here we 
just quote the relevant results for the calculation of the 
vacuum mean values of the EMT and current vector. 

The electric field acting during the time T creates a 
considerable number of pairs only in a finite region in 
the momentum space. One can introduce a cutoff K >> 
max{l,m 2 / \qE\) such that, for v / \qE\T/2 > K, one 
needs to consider only the region 



1/2 



Ip±I < V\qE\ VW\t-k 



-T/2 + K/^\q~E\ < Pl /qE < T/2 - K/y/\qE\ (12) 

in the momentum space. After the cutoff, physical quan- 
tities that are expressed in terms of K p will, in general, 
depend on K. But if T is big enough, so that 



^AqE\T > K > max {l, m 2 / \qE\) 



(13) 



the dependence on K can be ignored. We shall suppose 
that the stabilization condition (IT3"|) holds true. In this 
case the differential mean numbers of created pairs have 
the form 



(14) 



which is the same for the case of the constant uniform 
electric field. 

Taking into account cutoff (fT2|) . we find that the to- 
tal number H of pairs created is proportional to d — 1- 
dimensional spatial volume V^-i) and can be expressed 
through the total number density n cr of pairs created 
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during the interval T as follows: 



(27T) 



T+\qE\- 1,2 0{K) 

d/2 



o[d/2]-l 



where the quantity r cr is often called the pairs production 
rate. 

It should be noted that the stabilization condition is 
written for arbitrary external field E. However, in what 
follows we are interested in strong electric fields 



E > E c => m 2 / \qE\ < 1. 



(16) 



In what follows, when ever we speak of the strong field 
case we will have in mind the condition (I16[) . In these 
cases the stabilization condition is simplified, 



y/\qE\T » K » 1. 



(17) 



The stabilization condition reveals another important di- 
mensionless parameter r, 



(18) 



In the strong T-constant field under consideration the 
qualitative supposition that T is big enough is equiva- 
lent to the supposition that r> 1. In what follows, we 
calculate different mean values in an approximation that 
is related to large r; we call such an approximation large 
t- limit. To explain the meaning of such a limit, we con- 
sider possible structures of the mean values that appear 
in our further calculations. The most general structure 
of the mean values has the form 



F + Flnr + F%t 



For 2 . 



(19) 



The large r-limit for the mean value (F) means the lead- 
ing term approximation for (|19p. having in mind the hi- 
erarchy F <§C Flnr <C Fit -C F 2 t 2 . For example, the 
r.h.s. in Eq. (fT4)) represents the term Fq, while r cr in 
Eq. (fT5|) represents F\ in the general formula (fP9|) . 

It should be noted that below we encounter cases where 
the large r-limit is defined by t = \/\qE\At, with some 
At, such that r > K > 1. 

It should also be noted that for large T, the in- and 
owt-solutions in the interval t 6 Int take the form: 

+<p P ,s{t) = CD I/ _( 1+s ) / 2 (±(1 - i)Q , 
±V»p,.(*) = C£U_ (1 _ s)/2 (±(1 + , (20) 

with the normalization constant 

C = (2tt)-^- 1 ^ 2 \2qEl 1/2 exp(-7rA/8) . (21) 



It is supposed that the measurement is carried out 
at some time after switching off the electric field, i.e., 
decoherence occurs after the electric field is switched 
off. Of course, one can consider the case when deco- 
herence occurs earlier, for example, at an instant tdec, 
tin < tdec < t ou t ■ Let us suppose, for example, that the 
interval tdec — Un satisfies stabilization conditions simi- 
lar to l[I5j). y/\qE\ (tdec - ti n ) > K, Then the differential 
mean number of pairs created in the large tdec — tin limit, 
when 

1/2 



|P-L| < Vk^l V\QE\ (tdec - U n ) - K 



t m + K/^\q~E\<p 1 /qE<t d ec-K/^/\q~E\, (22) 

is given by (|14|) and the total number density of particles 
created is r cr (tdec — Un)- F° r the time of decoherence 
tdec, which is sufficiently close to t out , (t out - tdec) /T < 
1, the difference of this density from the density n cr can 
be ignored and the interpretation of particles at tdec as 
final out-particles already makes sense. 

However, if the interval t ou t —tdec is comparable with 
the interval tdec — U n then it is necessary to consider the 
further evolution of many-particle state given initially by 
the distribution (| 14[) at the time tdec, which is a com- 
pletely different task, see for example (25[. In this case, 
the differential mean numbers of additional pairs created 
by the external field during the interval from tdec to t out 
is given by Eq. (59) of the work [2|| as follows: 



AJVp = H p 



where 7V p +) (m) = iV p '(in) 
numbers of initial particles and antiparticles at tdec, given 
by Eq. flU]). If N { p +) (in 



1 - 



(AK+)(m) + AT p -)(m))] , (23) 



(-)/ 



Hp are differential mean 



A^p ] (in) 



> 1 then A7V p is 



negative, that is, the annihilation of the existing pairs of 
particles and antiparticles results from the electric field. 
Integrating and summing distribution (1231) over the quan- 
tum numbers under conditions (|22|) , we obtain the addi- 
tional number density of pairs 



An c 



1 



j2-d/2 



exp 



\qE\J 



f (tdec tin) 



(24) 

density of pairs created at tdec, 
the additional density An cr , and the 
number density of pairs created from vacuum states after 



Summing the 

^* (tdec ti n ) , 



tde 



(tout — tdec), we find the following total number 



density of pairs created during the interval T, 

n cr (tout - tdec | *dec - *in ) = n° r + An cr . 

Thus, if decoherence occurs at tdec, the final number den- 
sity of pairs at the time instant t ou t differs from n cr by 
An cr . In the case of strong electric field, t^W <C 1, 
An cr for d = 4, An cr < for d = 2, 3, and An cr > 
for d > 4. We see that the intermediate decoherence at 
the early time significantly reduces the measured den- 
sity of final particles in low-dimensional systems, and in- 
creases in high-dimensional systems. Decoherence can 
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occur once or many times during the interval T. The 
choice of a suitable model of decoherence depends on the 
physical nature of the phenomenon. 

It is useful to comment the case when the initial many- 
particle state is the thermodynamical equilibrium of non- 
interacting particles at the temperature O with the chem- 
ical potential jj, c h, 

N^(in) = {exp [(Lu m - n ch ) /0] + ly 1 , 

where uji n is given by Eq. (|9]). Then the differential mean 
number AN p of additional pairs created by the external 
field during the interval T, has the form (l23|) , where H 



is given by (fT4"l) . Taking into account the cutoff (|T2l . one 
can find the total number of additional particles created. 
In consequence, one can see that at low temperatures, 
(to — fi c h) /© 3> 1, such a total number differs from the 
zero-temperature result (|15p by a next-to-leading term 
that is not essential in the large r- limit, see [251 ] . We 
consider the system at high temperatures when all the 
energies of the created and accelerated particles are much 
lower than the temperature 0, hence \ qE \ T/0 <C 1. One 
can extract from results of the work [25[ , that the total 
number of additional particles created under such a con- 
dition is much less than that in (|15[) . However, we stress 
that at high temperatures, the polarization effect in the 
current and energy densities of the initial gas of charged 
particles given byEq. (127) in QjJ, is much stronger than 
effect from the pair creation given by Eq. (72) in [Toj . 
This could lead to screening the electric field before the 
effects of pair creation manifest themselves. One can see 
that the vacuum contributions dominate in comparison 
with contributions due to the low temperature and par- 
ticle density in the initial state. That is why we restrict 
ourselves to the case of the vacuum initial state. 



III. ENERGY-MOMENTUM TENSOR AND 
CURRENT VECTOR 

A. Observables and Green functions 

We consider here various singular functions of the 
Dirac field. With the help of these functions, different 
physical quantities can be calculated. For example, mean 
values © can be calculated with the help of the so-called 
«n-propagator: 

Sf n (x,x') = i{0,in\T^(x)^(x')\0,in) . 

In turn, this propagator can be determined via the in- 
solutions as follows 

Sl(x, x') = 6{t - t')Sr n {x, x') - 6{t' - t)S+(x, , 

Sf n (x,x') =i / ±%A X ) ±$p,r( X ') ■ ( 25 ) 

r 

Using relation ([2]) and properties of the g coefficients, one 
can divide the propagator (|25[) into a sum of two terms, 

Sf n (x, x') = S c (x, x') + S p {x, x') , (26) 



where the first term is the causal (Feynman) propagator, 

S c (x,x') = i(0, outlT^ix^ix'^OJn)^ 1 , 

c v = (0, oui|Q, in) . (27) 

Here c v is the vacuum-to-vacuum probability amplitude. 
This propagator can be represented as follows: 

S c (x, x') = 6(t - t') S~ (x, x') - 6(t' - t) S + (x, x') , 



S (x, x')=i y dp + 1p v ,r ( 



x )9{+\ 



S + (x,x') = i dp^2-ip p , r (x) g(-\ ) 1 ip P A x ')- 

J r 

(28) 

The second term in (|2"6"]l has the following form: 
S p (x,x')=i /dp^_Vw(£) 

J r 

x [ 5 ( + |-) 5 (_|-)- 1 ]* + %A X ') • (29) 

The current density, j£.(t), and EMT, T^(t), of cre- 
ated particles at t > t out , are expressed via the mean 
values of the normal form of j M and T^ v operators with 
respect to the out-vacuum. Namely, 

«) = 0'"(t)>-0'"(t)>»t. 

{f{t))out = (0,out[f\0,out), 

i;(t) = W))-(W)u. 

(T^(t)) out = (0,out\T^\0,out) . (30) 

These mean values can be calculated with the help of the 
so-called out-propagator 

S c out {x, x') = i(0, out\TV(x)V(x')\0, out) . 

Similarly to the case of the in-propagator, one can relate 
this propagator to S c (x, x 1 ) as follows: 



S c out (x, x') = S c (x, x') + Sf(x, x 1 ) , (31) 



where 



S p (x,x') = -i [dpJ2 + Vw(z) 

J r 

xg{ + \ + T l 9{ + \-)-%A x ')- (32) 

The quantities Q and (|3T)]) are real valued and can be 
represented as: 

(,f(t)) =Re(r(t)) c + Re(r(t))P, 

(T^(i)) = Re (T^(t)) c + Re (T^t))* , (33) 

(j^t)) out =Re(r(t)) c + Re(j^t)y, 

(T^(t)) out = Re (T^(t)) c + Re {T^{t)Y . (34) 
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where 



{f{t)Y>™ = iq tr [<fS°**(x,7>)} \ x=x , , 
(T,At)) c - p - p = * tr *')] \ x=x , . 

V = 1/4 [ 7 „ (p„ + i?) + 7 „ (p„ + p;;)] 



(35) 



Here tr stands for the trace in the 7-matrices indices and 
the limit x — > x' is understood as follows: 



tr[- • • (x,x')] x=x , = - 



lim tr[- • • (x, x')] 



+ lim tr[- • • (x, x')]] x =3( 

t-yt'+O 



The representations (|3"3")l and (jM)) imply that 

«)=Re<jW-Re<jW~, 
T£(t) = Re (T^(t))f - Re (T^{t)) p . 



(36) 



Note that the mean current (j M (i)) and the physical 
part of the mean value (T^it)} are zero for t < U n , when 
the electric field is zero. We are only interested in these 
mean values for large T and for t 6 Int, when the time t 
from the latter interval is sufficiently large, 



qE(t - U n ) > K » max{l,m 2 /|gP|} , (37) 



where K is the cutoff introduced before in Eq. (|22|) . Then 
it is sufficient to use the in- and out-solutions in asymp- 
totic form ([20]) for the functions S c ' p (x, x') defined above. 
At late times, t > t out , the solutions ± tp P:r (x) reduce to 
free plane waves in accordance with Eq. (J9)). 

Some components of (j^ (t)) c ' p ' p and (T^„ [t)) c - p - p are 
finite (do not have any T-divergences) as T — > 00. In the 
expressions for these components, we can use asymptotic 
(as T — > 00) forms of the singular functions S c,p ' p (x, x'), 
that is, the functions S c,p,p (x,x') in the constant electric 
field. For such singular functions, the so-called Fock 
Schwinger proper time representations hold true [26( | (it 
should be noted that the functions S p (x, x') and S p (x, x') 
defined above coincide with the functions —S a (x,x') and 
— S p (x,x') respectively, used in the article [26]). In 
the Fock-Schwinger representations the causal propaga- 
tor S c (x,x>) defined by Eq. flSS), S p [x,x') in and 
S p (x,x') in have the following integral representa- 



tions: 

S c - p - p {x, x') =( 7 P + m)A c ' p ' p (x, x') , 
dsf(x,x',s) , 



A c {x,x') 
A p (x,x') 

A p (x,x') 



ds f(x, x 1 , s) 



Q(xi-x[) dsf(x,x',s), 
Jr 3 +r 2 -r p 

ds f(x, x , s) 



Q(x[ - X!) / dsf(x,x',s) , 

Jr 3 +r 2 -r p 

(38) 



where the Fock-Schwinger kernel f(x, x', s) reads: 
f(x, x',s)= exp (-i^F^s) /W (x, x',s), 

r '{x,x ,s) 



{4iri) d / 2 smh(qEs) 
1 



e iqA e~ im s 



x exp 



4/ 



(a; — a;')gPcoth(gPs)(a: — a;') 



(39) 



coth(gPs) is the matrix with the components 
[coth(qPs)]^„, and A = (x + x' Q )(xi - x' x )E/2. 
All integration contours in the s-complex plane are 
shown in Fig. [TJ The integral along the contour r c , that 
is, along the real positive semiaxis, corresponds to the 
well-known Schwinger's representation of the Feynman 
propagator. 

Im s 



O 



Re s 



■2\qE\ 



\lE\ 



FIG. 1. Contours of integration r2,r3,r c ,r p . 
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The only singular points of the kernel f(x,x',s) in 
the lower half-plane outside of the origin are s n — 
—mnj \qE\, n = 1, 2, . . .. 



B. Peculiarities in odd dimensions 

We shall see that in odd dimensions there exist current 
components related to the so-called Chern-Simons term 
of the effective action. The components (j 1 (t)) p ' p are 
ill defined in the proper time representation due to the 
divergence in the limit T — > oo. They will be treated in 
detail at finite T in Subsection IIIIEI All other current 
components (J** (t)) c,p,p are finite as T — > oo and can be 
determined with the help of Eqs. (J3SJ. We first consider 
the case of d = 3 dimensions. In this case, there are two 
non-equivalent representations for the 7-matrices, 
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a 3 , 7 1 



r 



= -i(±l)a\ 



(40) 



where a % are the Pauli matrices and the signs ±1 cor- 
respond to different fermion species, which we call as 
±-fermions, respectively. By using the formula 

cxp (-i^cr^F^ = cosh(qEs) + 7V siah(qEs), 

one finds that 

tr [^mf{x, x, s)] = -i(±I)2msinh(g£ , s)/ (0) (x,x,s)(^2- 

Inserting this expression and f^\x,x',s) given by 
Eq. (|39"|) into Eqs. (|38|) and taking into account that the 
integral along the contour + T2 
tribute to quantities given by Eqs. 
following result: 



T p does not con- 
SJ, we obtain the 



e 2 E 



-in 



[J[ )) T 4n*/* \2' \qE\J ' 



(41) 



where r(l/2,x) is the incomplete gamma function. In 
d =/= 3 dimensions we find that 



(42) 



We see that in d — 3 dimensions there are nonzero com- 
ponents (f{t)) c f' p = 5^ 2 {j 2 {t)) c ^ p that are orthogo- 
nal to the electric field direction. The nonzero value of 
(j^(t)) c indicates that the Chern-Simons term is present 
in a properly regularized effective action with odd num- 
ber of fermion species, see [13, [28j]. The signs "±" in 
Eq. (|4T|) are opposite for each of the two possible fermion 
species. The absence of similar components in expres- 
sions (|4*2"j) in higher (d > 3) odd dimensions is related to 
the fact that we consider a special case where the mag- 
netic field is absent. Note that in arbitrary constant field 
such current components are orthogonal to the electric 
and magnetic field, {j^{t)) c f' p F^ = 0, and proportional 



to the product of all eigenvalues of the field tensor F MV , 
see the appropriate expression for tr [j^rnf (x, x, s)] in 
26]. Separating the transverse components, we repre- 
sent the vacuum current density in the final form 



<j"(f)>x = (f{t))[ 



Mf(t))l 
+ {f{t))i, 



(43) 



where finite quantities (j M (i))^ p are given by Eqs. (|4Tj) 
and (fi2]l. whereas the component (j M (i))j| = £ M ,i(j 1 (*)) P 
directed along the electric field must be treated for finite 
T. One can see from Eqs. ([36)) and (j4T|) that the trans- 
verse current of created particles is absent, j 2 r (t) = 0, 
and this conclusion holds true for arbitrary constant field 
in higher odd dimensions. Then the term (j l *(t))± repre- 
sents transverse vacuum-polarization current, in partic- 
ular for d = 3 it has the form 



±5 



'fj.,2- 



e 2 E ( 1 7rm 2 \ 



47r 3 / 27 I 2' \qE\ ) ' 



(44) 



where 7 (1/2, x) is the incomplete gamma function. 

The expression for j_ that can be extracted from 

(l4"TT) has previously been obtained in [2?], while the 
expression for (j fl (t)) p ^ p that can be extracted from (|4Tj) 
is new. It should be noted that the expression for the 
nonequilibrium Hall conductivity at long times given by 
(gH is close to the one given by Eq. (7) in [2{|. The lat- 
ter was obtained in the framework of one-particle WKB 
calculations, using analogy with the Landau-Zener tun- 
neling. Note that it is quite difficult in the framework of 
the one-particle theory to distinguish between the cur- 
rent of real particles, which remains after switching off of 
the electric field, and the pure vacuum-polarization cur- 
rent, which disappears in this case, although there is no 
dissipation in the model. We stress that the latter cur- 
rent disappears because of the switching off the electric 
field. Thus, we see that using in this problem analogy 
with the tunneling, or with the pair-production rate by 
Schwinger, has to be justified. In fact, our exact result 
(|4"4"|) . obtained in the framework of QED, is such a justi- 
fication of the one-particle calculation [29(, at the same 
time it shows the limits of its validity, namely, it holds 
true only before the switching off the electric field. The 
quantity (j 7i (i)) p L ' P is exponentially small for weak electric 
fields, m 2 / \qE\ > 1, in which case (j 2 {t)) w (j 2 (t)) c . On 
the other hand, in the strong-field limit, m 2 
both contributions are comparable, (j 2 {t)) p 
and then 



^ 2 /\qE\ « 1, 



e 3 / 2 Em 



(45) 



The conductivity that follows from eq. (j4"5"j) coincides 
with the one given by Eq. (9) in [29]. 

As it is known, the term (j 2 {t)) c is related to the stan- 
dard effective action and therefore to probability ampli- 
tudes of processes, while the term (j 2 (t)) presents a con- 
tribution to mean values which, in general, are quite dif- 
ferent from such amplitudes. We see that (j 2 (t)) c and 
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the probability amplitudes remain unchanged as m — > 0. 
In contrast to the behavior of quantity (j 2 (t)) c , the mean 
value (j 2 (t)) tends to as m —> 0. It is important to note 
that these constant values of both (j M (i))5_ an d 
are obtained in the limit T — > oo, that is, Eqs. (|4"Tj) 
hold true for instants of time t before the electric field 
is switched off. At early, t < U ni and late times, t > tout, 
the T-constant electric field is absent, and one can see 
from the exact formula (|35j) that the vacuum polarization 
current vanishes in this case, (j t _ t (t)) c — (j M (i))^ = 0. 
These results can be generalized to the case of massless 
fermions in an electric-like constant electromagnetic field 
in higher odd dimensions, when the electric field cannot 
be removed by the Lorentz transformation and all eigen- 
values of the field tensor F av are different from zero. In 
this case when t < tout, the transverse component of vac- 
uum polarization current (j fJ '(t)}j_ is nonzero, however, 
the total transverse mean value is equal to zero. 

Thus, the vacuum polarization in the electric field in odd 
dimensions is qualitatively different for mean values and 
amplitudes of processes. This is what distinguishes sub- 
stantially the case of electric-like field from the case of 
magnetic-like field (in the latter case the electric field 
can be removed by a Lorentz transformation). 



C. EMT 

Using representations ([55)1 and (|3"5)l . we obtain compo- 
nents of the EMT as follows: 



(T M!/ (t)) c ^ = 0, p + v, 



(46) 



Re(T 00 (t)) c = -Re<T u (*)) c = E dRe ^ _ ReC (t) , 



dE 

Re{T tl (t)) c = Re£(t), i = 2,3,...,D, 



(47) 



where 



„ 1 f 00 ds ,, 

£ = - / —trf{x,x,s) , 

trf(x, x, s) = 2 [d/2] cosh(qTs)/ (0) (a;, x, s). 



(48) 



It should be noted that the diagonal elements {T uu (t)) p 
are ill defined in the proper time representation due to 
the divergence at the limit T — > oo and will be treated 
at finite T in Subsec. IfflEl The quantity C in (gSJ can 
be identified as the non-renormalized one-loop effective 
Euler-Heisenberg Lagrangian of the Dirac field in an uni- 
form electric field. It is the density of the one-loop effec- 
tive action W, 



W = / Cdtdx = 



(49) 



that is defined in general via the vacuum-to- vacuum am- 
plitude (|27p. Its imaginary part represents the vacuum- 
to- vacuum probability, as follows: 



It is a global physical quantity. It is free of ultraviolet 
divergences, because creation of pairs with infinitely large 
momenta is suppressed. The Bogolyubov transformation 
between in- and out- creation and annihilation operators 
allows one to relate c v with the g coefficients from ([2]) 
QUI: 

c v = exp {trln g(_|~)*} . 

For the T-constant field, when T satisfies the stabi- 
lization conditions (fT5|) . the probability (150)) can be ex- 
pressed via total number of particles created (fl~5|) in the 
following form: 



oo 



ln|c,,| 2 = -pK, p=J2(n+ir d/2 exp( 

n=0 



2 /\qE\) , 

(51) 



-2ImW 



(50) 



see [24|. It follows from dHIJ) that H ~ TV {d -i), then 
ln|c u | 2 is proportional to the spatial volume V(d~i) and 
the field duration T. It can be verified that the result 
(f5Tj) can be obtained with the help of Im C (|4"8|) , provided 
that the integral over the time t is identified with the field 
duration T, Jdt = T, Thus, Im W is finite for finite 
values of V(d-i) and T. Note that in the strong-field case, 
m 2 / \qE\ <§; 1, for d > 3, the leading term in the r.h.s. of 
(|51[) is the ordinary Riemann zeta function, p — Qr (d/2). 
For d = 2, the leading term is p = In (\qE\ /m 2 ). That 
means that in d = 2 dimensions the vacuum-to-vacuum 
probability of massless fermions in quasiconstant electric 
field is ill defined. In this case, the mass term m ^ has 
to be considered as an infrared cutoff. 

Both the vacuum mean values of the total current den- 
sity and of the EMT are represented as sums of two con- 
tributions: Re(j^(i)) c (resp. Re(T M „(i)) c and Re£(i)) 
that can be associated with the vacuum polarization, 
while (j^(t))P =JZ r {t) + (?"(*)>* (resp. (T^(t))P = 
T^it) + (T /1 v(t)) p ) can be associated with the pair cre- 
ation due to the term j£.(t) (resp. T^(t)). The latter 
follows from the fact that S p (x,x') with all its deriva- 
tives, given by Eq. (|2"9"|) , is exponentially small due to 
the smallness of |<?(+|~)|, when the electric field is weak, 
m 2 1 \qE\ I ^> 1. Of course, in general, such unambigu- 
ous division of physical quantities due to particle creation 
and the vacuum polarization is not possible. However, it 
can be done in some specific cases and for some specific 
quantities. It is clear that j£ r (t) and T^{t) depend on 
the history of the process and retain their latest values 
at t > t out ■ On the other hand, we see from Eqs. (|4"Tj) 
and (|4T|) that the nonzero real parts of the quantities 
(^(i)) c , Re(T u „(t)) c , and are time-independent 

for the time t € Int. At early t < ti n and late times 
t > t out , we have (j /i (t)) c = (j^it)) 1 ^ = and the mean 
{Tuv{t)) c is reduced to its free value for E = 0. Thus, we 
see that (j^{t)) c , W(t)) p ± , and Re(T M „(i)) c depend on 
the electric field at the time t, but do not depend on the 
history of the process, that is, they are local quantities 
and represent the vacuum polarization contribution. We 
continue to study each of these local and nonlocal terms 
independently in the following sections. 
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D. Mean values: renormalization 



The integrals (|47|) are divergent due to the real part 
of the effective Lagrangian (|48p which is ill defined. This 
real part must be regularized and renormalized. In low 
dimensions, d < 4, Re£ can be regularized in the proper- 
time representation and renormalized by the Schwinger 
renormalizations of the charge and the electromagnetic 
field In higher dimensions, a different approach is 

required. Here one can give a precise meaning and cal- 
culate the one-loop effective action using zeta-function 
regularization [3l[ . An application of this method to the 
case of a uniform magnetic field and self-dual field in ar- 
bitrary dimensions is described in detail in [32] . It can be 
shown that for d < 4 such a renormalization is equivalent 
to the above mentioned Schwinger's renormalization, see 
32] . Let us consider the application of this technique for 
the case of a constant uniform electric field in arbitrary 
dimensions of interest here. 

First, we remind that the effective action W from ([49]) 
can be represented as the functional determinant, 



W = — IndetM 2 , M 2 
2 



ft is ■ 

(52) 

The operator M 2 becomes elliptic when it is continued 
to the Euclidean space, M 2 —> M 2 , by means of the 
replacements t — > —ir), do — > id^, and qE — > iB. Then 
the functional determinant is well-defined with the help 
of the zeta-function Q d ' (s) in d dimensions, 



IndetM 2 



dC {d) (s) 



ds 



C W (s)=T- 1 ( S ) 



s=0 

duu s ~ x 



K(u) 



(53) 



where K (u) is the heat kernel, 

K(u) = J drjdxtrf E uci(x,x,u) 



fEucl(x,X,u) 



exp 



(-«m7c 2 )|i 1 x) 



Here the quantity [i is a renormalization scale, which is 
introduced to keep the zeta-function dimensionless. The 
dependence on /i of the functional determinant corre- 
sponds to a finite renormalization. 

The effective action W can be written as the fol- 
lowing integral over the Euclidean space time volume: 
W = —i J Cdrjdx. It is real when it is continued to the 
imaginary electric field qE — > iB, so that W = W\ qE=iB , 

ImW = 0. Note that for d > 3, W coincides with the 
effective action of the magnetic field B/q. Consequently, 
we obtain a regularized and finite form of the effective 
Lagrangian of a constant electric field in an arbitrary d- 
dimensions as follows: 



Re£ 



reg 



Re £ 



B=-iqE 



C = -2 



1 dC (d) (s) 



ds 



s=0 



ft" 1 

"(d) 



where 0(d) = f d-qdx is the Euclidean space-time volume 
in d dimensions. 

An explicit form of the quantity txf Eu ci{x, x, u) can be 
extracted from the quantity tr/(x,x, s) (given by (|39[) 
and (1481)) as follows: 



trf Euc i(x,x,u) = - txf 



ill 



(55) 



qE=iB 



Then we have: 



c (2) w 



n 



(2) 



i 



2tt F(s) 



ft 



(2) 
K2) 



4tt J 



^coth 



e-^/^du. 



x(W (s-d/2 + 1) 



r(s-d/2+l) 
for d > 2 . 



(56) 



If B > 0, the function Q 2 '(s) can be written in terms 
of Hurwitz zeta-function as follows [3ll 13211 : 



C (2) (*) = 



^(2) B 
2-r 



^(2)f 



■IB 



■m 

for m 7^ 

The Hurwitz zeta-function is defined as an analytic 
continuation of the series 



Cr{s) , for m — . 



Ch(s;x) 



n=0 



(n + x)~ 



Re s > 1 , 



to the entire complex plane of s. One can see that 
Cr( s ) = Ch(s, 1) is the ordinary Riemann zeta-function. 
For odd d, and taking into account that r~ 1 (s) m s at 
s w 0, we get from ((5E)) that 



dC (d) (s) 



ds 



= r(,)c (d) (*) 



(58) 



By using Eq. ([54]) we find a finite expression for Re C re g 
in arbitrary dimensions. Then we have to implement 
an additional finite renormalization of the cosmological 
constant and the electric charge [32j. For d > 4, a fi- 
nite renormalization of some high-dimensional quantities 
could be needed. Thus, we obtain the final form of the 
renormalized effective Lagrangian Re C ren . The corre- 
sponding final forms for C and for the B field in d = 2, 3, 4 
dimensions are treated in detail in [32l |. 

We are interested in the case of a very strong field, 
m 2 / \qE\ « 1. In this case the leading contribution to 
Re£ re „ is given by Re£ reg in |54|) . where m — > 0. Using 
([5o]) . (|57|) . and ([55]) . we find that this contribution has 



the form 



In 



C (d) (Q) 



r( s ) 



C (d) (*) 



B=-iqE 



B=-iqE 



(54) In particular, 



for even d, 
, for odd d. 
(59) 
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Re£ r 



Re 
Re 



2^(1) Cfl(|) 



_ J?£i 

— 4 



Re 



-^rln 



B=-iqE 



B=-iqE 

1 liT 



igg| 3/ V , 



for d = 2 , 
for d = 3 . 
for d = 4 . 



(60) 



For d = 4, the obtained Re£ ren coincides with the al- 
ready known result @. For the magnetic field, the lead- 
ing term in the form of the second Eq. (1601) has previously 
been obtained in [28[ for d = 3. Note that despite the 
fact that the nonzero vacuum current (j M (i)) c is present 
in d = 3 dimensions, the Chern-Simons term vanishes for 
the constant field strength [28J. In general, we have for 
a very strong electric field that 



Re C r 



\qE\ d l 2 , d^4n, 
\qE\ d/2 \n(\qE\/fi 2 ) , d = An . 



In contrast to the electric field case, where the logarith- 
mic factor In (\qE\/fi 2 ) appears only for d — 4n, in a 
strong magnetic field it is present in any even dimension: 
Re£ re „ - B d l 2 MB l[i 2 ) for even d, and Re£ re „ - B d l 2 
for odd d. In the framework of the on-shell renormaliza- 
tion of massive theory, we have to set fj, = m. 

Thus, we have obtained the renormalized mean values 
of EMT components in the following form: 



Re(T 00 (t)) c ren = -Re(T n (t)) c re1 
,9Re£ re „(t) 



= E 



ne(T«(*))J 



dE 

R,G £iren (^) ? 



R,6 £ reTl (£) , 

= 2,3,..., A (61) 



where C (t) ren at t € Int is given by (|50[) . Thus, in the 
strong-field case, the quantities (|6"Tj) have the following 
behavior: 



Re(T^(t)) c re 



\qEf 2 In {\qE\ / ^) , d = 4n . 
|g£| d/2 , d^4n, neN. 



(62) 

At early t < U n and late t > t otlt times, we have 
ReC(t) ren — and Re(T M1/ (£))£ en = 0. Taking into ac- 
count (|6ip. we find the final form for the vacuum mean 
values of the EMT 

(T^(t)) ren = Re (T^(t)) c ren + Re (T IMU (t)) p , (63) 

where, according to (|4"6")l , the off-diagonal elements 
are equal to zero, and the diagonal elements of 
Re (Ttj,u(t))ren are gi yen by dH]); the diagonal compo- 
nents of Re (T^{t)) p are studied in detail in the next 
section. 



E. Mean values: pair-creation contributions 

We know that the mean values (j 1 (t)) p,p and 
{Tf i f l {t)) p,p are divergent as T — > oo. Therefore, in the 



representation ([331) . they have to be considered always 
at finite T . Let us evaluate these quantities in the large 
r- limit approximation with r = \J\qE\(t — U n ). Consider 
the time t G Int for which the interval t — ti„ satisfies 
condition Q37p and the cutoff approximations 



P±|<Vp1 y/\qE\(-t-tin)-K 



1/2 



U n + K/y/\qE\ < Pi/qE < t - Kj \J\qE\ 



(64) 



hold true. In Eqs. (164)) it is taken into account that phys- 
ical observables in the time moment t € Int are affected 
by the electric field that acted for the time [£;„ , t] . Let us 
call the corresponding region in the momentum space by 
f2(i). Using the transformation ([2J, one can represent the 
propagator S p in Eq. (f2T)f in terms of owi-solutions. In 
the asymptotic regime, as z — > oo, the Weber functions 
have the following asymptotic expansion: 



D v (z) ~ z^exp l 



74) 



' N 

E 

,n=0 



( 2)n (2 



n\(-z 2 /2) n 



which is valid for |argz| < 37r/4. Then, keeping only 
the zeroth (n = 0) term in last equation, we obtain for 



S p (x,x r ) =( 7 P + M)A p (a;,a;'), 
A p (x,x r ) =-i J dp2\qEt-p 1 \X p e 
x [ + ^ P ,- 1 (i) + ^ ! _ 1 (f) 

+-<P P ,+i(t)-<p* t+1 (f)] 



?'p- (x— x ) 



(65) 



(66) 



Considering the large r-limit in representations (|35|) . the 
domain of integration in (|66p can be restricted to the 
region f2(t) described by the inequalities (|64|). Using 
Eq. (fl4|) for the differential mean values H p , we obtain: 



A p (x,:r/) 



h \\( x h x '\\) = ~j~i c 

\qE\^ 

x exp 



K/y/\qE 
,<Pl(*i-«l) cos 



^j.(xj.,Xj_)ft||(i||,a;j|)dt, 
1 



[€(tO a -£(*) a 



1^1 



(x ± -xl) 2 | gj E| 

47T 



(67) 
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where p\ = qEt. 

It follows from Eqs. {35) and {57) that 

(j\t)Y = -iq2VM P x K>{x,x')\ x=x , , 
(T„(t))P = l 2W^P^(x,x')\ x=xl . (68) 

Integrating over p\ , we obtain the following result in the 
large r-limit (with t = ^J\qE\(t — ti n )): 



= 2e sgn (E) r cr It - t m + \qE\-^ 2 (K) 



(69) 



(T Q0 (t)) p = (TnW) p 

- |<?£| r cr [t-t m + \qE\- 1/2 0(K)\ 2 , 
(Tn(t))* = 7r-V cr {in [vp(< - i in )] + O (InK)] 



for i = 2, 3, . . . , D, 



(70) 



where r cr is given by (|15). We see that all the leading 
contributions given by {69) and (|70) are real. The quan- 
tities (j 1 (t)) p and (T atl (t)) p depend on the time interval 
(f — Un) of the electric field action, showing that they are 
global quantities. 

Now we estimate, for t w t out , the current density 
and EMT of created particles, jg r (t) , T^(t), given by 
Eqs. ([55) . At this time instant, the solutions ±- 0p,r(^) 
reduce to free particle plane waves in agreement with 
Eq. {9). Thus, taking into account representation {32), 
one can see that the quantities and (T U u(t)) p 

from (|55) can be neglected in the large r-limit. Using 
Eqs. {69) and {70), we obtain: 



(71) 



where Eqs. {3D, (|42l) . and {35) were taken into account. 
For t > t out (after the electric field is switched off), the 
quantities j£ r (t) and T^(t) are constant and retain their 
values at t out ; for this reason, we always have to set t—ti n 
= T in such cases. The renormalized vacuum polariza- 
tion contributions in expressions {43) and {63) vanish in 
the absence of external field. Therefore, for t > t ou t, the 
vacuum mean values (j^it)) and T av (t)) r en represent the 
mean current density and EMT of pairs created by the 
complete T-constant electric field: 

(j" (*)> = far (tout ) , (T^{t)) ren =T^(t out ), t> t out . 

(72) 

We see that the current density of created particles 
jg r (t) is directed along the direction of the electric field. 
In d = 3 dimensions, in contrast to this almost obvious 
property, for t in < t < t out , the mean current density of 
each massive fermion specie deviates from the direction 
of the electric field. Indeed, 



0' 2 (*)> 



1 nm 2 \ exp (nm 2 / \qE\) 



(f(t)) ± 2 7 U \qE\) yW\(t 



(73) 



where yy\qE\Tj -C 1 according to stabilization con- 
dition (|13|) . One finds a similar deviation of the mean cur- 
rent density of each massive fermion species in higher odd 
dimensions for t in < t < t out in the case of an electric-like 
constant electromagnetic field, when all eigenvalues of the 
field tensor F av are different from zero. The total mean 
current density of an even number of fermion species is 
directed along the direction of the electric field, since the 
contributions of the ±-fermions differ only in sign. 

One can see from Eq. ([55) that r cr T is the total number 
density of pairs created and accelerated during the time 
T to velocities nearly the speed of light. It coincides with 
the quantity n cr obtained in a different manner in {15). 
The quantity T§q — \qE\ Tn cr is the mean energy den- 
sity of pairs created at any time instant t £ Int with zero 
longitudinal kinetic momentum and then uniformly ac- 
celerated to kinetic momenta from zero to the maximum 
\qE\ T, so that \qE\ T/2 is the mean kinetic momentum 
per particle. The energy density Tqq is equal to the pres- 
sure T^l along the direction of the electric field. This 
equality is a natural equation of state for noninteracting 
particles accelerated by an electric field to relativistic ve- 
locities. The momentum density of created pairs T§[ is 
zero due to the symmetry between particle and antipar- 
ticle distributions. 

The vacuum mean values {33) and ([53) for t 6 Int 
are sources in equations of motion for mean electromag- 
netic and metric fields respectively. It should be noted 
that only when the time t is sufficiently close to tout, 
[tout — t) /T <C 1 , the differences between the densities 
(j 1 (t)) p and (T l j, l/ (t)) p and the respective densities j\ r (t) 
and Tj%, (t) of final pairs created, can be neglected and 
the interpretation of particles at t as final owt-particles is 
correct. 

In the general case when the time t is not close to 
t ou t, there is an essential difference between the defini- 
tion of the vacuum at t < t out and the final vacuum 
state \0,out) at t out . That is why the quantities (j 1 (t)) p 
and (T fl „(t)} p have nothing to do with characteristics of 
final owf-particles. They present contributions to mean 
values due to vacuum instability which depend on the 
history of the process, that is, they are global quanti- 
ties, in contrast of the local quantities (j At (£)) c , (j^(f)Y\_, 
and Re (T M y(t)) c . In the general case, r cr in Eqs. ([55) 
and (|70p is the total number density of excited states per 
unit of time. For example, the longitudinal component 
of the mean current density (j 1 {t)} p increases linearly as 
(t—ti n ) grows, since the decoherence does not take place 
for t < t out . We note that (j'^W )) « jl r (t out ). How- 
ever, it maybe not so if the decoherence starts earlier, 
for example, at the time instant tdec, td ec < t out , and 
tout — tdec is macroscopic; see discussion in the end of 
Section [III In this case, the quantity (j 1 (t ou t)} decreases 
significantly in low-dimensional systems, and increases in 
high-dimensional systems. 

We can compare contributions from the vacuum insta- 
bility with contributions from the vacuum polarization. 
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Of course, all contributions due to pair creation in ex- 
pressions (jf?3")) are exponentially small for the weak elec- 
tric field, m 2 / \qE\ >• 1, so that the vacuum polarization 
terms are principal. We are interested in the strong-field 
limit, m 2 / \qE\ <C 1. In such a limit, we obtain from (|70|) 
that 

(Too (t)) p = (T n (t)^ ~ \qE\ d/2 \qE\ (t - t m f , 

(T u (t)f ~ \qE\ d,2 \n[^/\q~E\{t-t in )\ , % = 2, 3, . . . , D, 

(74) 



where the large dimensionless parameter \f\qE\{t — t{ n ) 
satisfies the stabilization condition (|37|) . Comparing 
the evaluation of the EMT components from (|74p and 
(|62p. we see immediately that, when d is not a multi- 
ple of four, the energy density of vacuum polarization, 
Re (Too (i)) rem i s negligible compared to the energy den- 
sity due to pair creation, (Too due to inequality 

If d = 4, the ratio \(T 00 {t)) P /Re (T 00 (t)) c | in a mas- 
sive theory with on-shell renormalization, u — m, is of 
the order \qE\ (t — t in ) 2 / In (\qE\ /m 2 ). In order to es- 
timate the allowed values of the logarithm in the latter 
equation, we have to have a physical model that describes 
the origin of the external classical quasiconstant electric 
field. In problems of high-energy physics it is usually as- 
sumed that just from the beginning there exists an uni- 
form classical electric field having a given energy den- 
sity. The system of fermions interacting with this field 
is closed, that is, the total energy of the system is con- 
served. Under such an assumption, we take into account 
that quantum electrodynamics in d = 4 dimensions with 
the strong T-constant external electric field can be con- 
sidered as a consistent model only if the backreaction 
due to pair creation is relatively small with respect to 
the background, which implies the following restriction 
from above: 



\qE\ (i-t m ) 2 « — 
2a 



(75) 



where a is the fine structure constant, see [10 . 
[33|]. We consider macroscopic time intervals, such 
that m 2 (t — ti n ) 3> I, then it follows from (f75|) 
that \qE\ /m 2 <C n 2 / (2a). Whence we obtain 
ln(\qE\/m 2 ) -C 6.5. On the other hand, accord- 
ing to condition (|75p . the maximum value allowed for 
\qE\ (t — ti n ) is two orders of magnitude greater than 
the restriction obtained for the logarithm. Thus, we 
see that in this case the quantity Re (T 00 (t))^. en is neg- 
ligible in comparison with (Too {t)) p - This result can 
be generalized to d — An dimensions. Note that 
the ratio of the transverse components of the pressure 

\(T u (t))P / Re (T u {t)) c \ is of the order In \y/\q~E\ (t - t i7 



if d ^ An and In y/\qE\(t - t m ) /In (\qE\ /m 2 ) if d = 

An. One can see, for example, in d — A dimensions, that 
these logarithms cannot be considered, in general, as re- 
ally big quantities due to restriction from above (|75[) and 



this evaluation can be generalized to arbitrary dimen- 
sions. Then, in general, none of these transverse compo- 
nents can be neglected. The evaluations may be differ- 
ent when another model for the external classical quasi- 
constant field is considered, for example, when there is 
an external source that supports a given external field 
strength. 



IV. MEAN CURRENT AND EMT IN 
GRAPHENE 

A. One-loop results in a given external field 

It is known that at certain conditions electronic exci- 
tations in graphene monolayer (just graphene, in what 
follows) behave as relativistic Dirac massless fermions in 
2+1 dimensions. The so-called Dirac model for elec- 
tronic excitations in graphene was developed first by Se- 
menoff in [l4l |. exploring results obtained decades earlier 
in the study of the conducti vity of graphite [Hj], see de- 
tails in recent reviews [HI, [T^- It was found that, at 
zero temperature and chemical potential (i.e., at the so- 
called charge neutrality point), low-energy electronic ex- 
citations in graphene are described in a tight-binding ap- 
proximation by the Dirac equation for massless particle 
in 2 + 1 dimensions, with the Fermi velocity vp ~ 10 6 m/s 
playing the role of the speed of light in relativistic par- 
ticle dynamics. In this section, we are going to explore 
such a correspondence, and consider applications of the 
results obtained in the study of the quantized Dirac field 
in an external background presented in the previous sec- 
tions to some problems of graphene physics that can be 
studied within the Dirac model. In fact, we are going to 
study the electronic transport in graphene at low carrier 
density and low temperatures when quantum interference 
effects are important. Accordingly, we shall restrict from 
now on to the massless case in 2 + 1 dimensions. 

First some comments about the Dirac model of 
graphene. There are actually two species of fermions in 
this model, corresponding to excitations about the two 
distinct Dirac points in the Brillouin zone of graphene, 
i.e., each species belongs to a distinct valley. The al- 
gebra of 7-matrices has two inequivalent representations 
in (2 + l)-dimensions, as described in (|4T)|). and a dis- 
tinct (pseudo spin) representation is associated with each 
Dirac point. There is no parity anomaly in the Dirac 
model, in particular, the sum of current densities {j^{t)) c 
(vacuum current contributions to the probability ampli- 
tudes for processes with photons) for the two fermion 
species, given by Eq. (l4"TT) . is zero. Note that the mean 
value of transverse vacuum polarization current, given 
by Eq. (|4"4l -. is equal to zero for each massless fermion 
species. For all other integral quantities, since interval- 
ley scattering can be neglected, the presence of two val- 
leys is taken into account simply by multiplying by the 
degeneracy factor 2. 

Furthermore, there also is a spin degeneracy factor. 
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The derivation of the Dirac model starts from a nonrela- 
tivistic Schrodinger equation for the conduction electrons 
in graphene, leading to another doubling of fields due to 
the (real) spin of the electron. As a result, there are four 
species of fcrmions in the Dirac model corresponding to 
graphene. The mean values that we have obtained in 
the previous section hold true for each of the Dirac fields 
independently. In order to find the corresponding mean 
values in graphene, one should first add the contributions 
from the two valleys, as discussed in the previous para- 
graph, and after that multiply by the spin degeneracy 
factor two. 

We consider an infinite flat graphene sample on which a 
uniform electric field is applied, directed along the plane 
of the sample. We assume that the applied field is the 
T-constant electric field studied in the previous sections: 
the field is suddenly switched on at some time tj„, act- 
ing then for a time-interval T, during which electron-hole 
pairs are created. We consider the case of zero tempera- 
ture and chemical potential, so that the Dirac model can 
be used, and an initial state with neither electrons nor 
holes. 

Under these circumstances, the Eqs. (fl5)) . (|43)l . and 
([55]) . multiplied by a degeneracy factor of four, describe, 
respectively: n g r ', the total number density of electron- 
hole pairs created by the electric field; (j 1 (t)} g , the mean 
longitudinal current density; and (T uu (t)) gi the mean 
EMT in the graphene. For the sake of comparison with 
known experimental results, we are going to use SI units 
and restore the Planck constant h in this part of the work. 
We get the following results: 

n cr = r cr T f r cr = ^-2 fatfyV* \ eE f 2 ; (76) 

(j 1 (t)) g = sgn(E)AAt, A = 2ev F r c ;; (77) 
(T 00 (t)) g = (T 11 (t)) g = e\E\v F r< g r (Atf, 

(T 22 (t)) g = (T 22 (t)r g + (T 22 (t)y g , 

(T 22 (t))c = fcr-C*(3/2)/2, 

(T 22 (t)) p g = hr^n- 1 In (Ve \E\ v F /hAt) , (78) 

where At = t — U n . These results hold true for all t that 
satisfy the stabilization condition (|37l) . which has now 
the form: 

y/e\E\v F /HAt > 1 . (79) 

There appears a time scale specific to graphene, 

At st = (e \E\ v F /hy 1/2 , (80) 

which plays the role of the stabilization time in the 
case under consideration. The vacuum polarization con- 
tribution to the mean value (Too(t)) g — (Tu(t)) g in 
(|78| is small due to Eq. ((79)) and is neglected, that is, 
(T 00 (t)) g = (T 00 (t))P and (T u {t)) a = (T u (t))f. 

At t ss tout, we have At « T = t out —tin and relations 
(I7T1) show that the mean values (j 1 {t)) g and (T a u(t)) p g 



hardly differ from the current density j\ r (t out ) and the 
quantity T^ r n (t out ) caused by created particles. In the 
general case, the quantity r g r is the number density of 
pairs of positive and negative charged states excited due 
to the constant electric field per unit of time. Only 
at i « t out , it can be treated as the production rate 
of electron-hole pairs. In the presence of a mass gap 
Ae = mv F , the rate r c g r is attenuated by a factor of 
exp [— 7r(Ae) 2 /e \E\ v F h] according to Eq. (TT5)) . In this 
case, the stabilization condition has general form (|57|) 
and the strong field condition reads (Ae) 2 /e \E\ v F H <C 1. 

The formula for the current density of created particles 
jcr(tout) that follows from expression (l77t agrees with 
the result obtained from the WKB approach [3, [l5|, [l9[ 
(the numerical analysis of (l9| in the frame of a non- 
equilibrium Green function (NEGF) approach, which is 
said to be a non-perturbative quantum mechanical ap- 
proach, is consistent with the semiclassical result). In 
condensed-matter physics such a method is known as the 
Landau-Zener approach (Note that the WKB approach 
is valid when the differential mean number H p , given by 
general expression (fTTj) . is small. However, the WKB 
approximation for H p coincides with the exact expres- 
sion (|14[) in the limit T — > oo). The time-dependence of 
mean current density (j (t)) g given by Eq. (1771) is con- 
sistent with the numerical solution of the first-quantized 
tight-binding model equations obtained for the ballistic 
case in the time interval At st < At < Ats, see Eq. (82) 
in Ref. [13]; there a factor 3 3 / 4 2~ 7 / 2 is replacing our 
2/7r 2 , but these factors are both equal to 0.20, numeri- 
cally. The AtB — 2nh(e \E\ a) -1 is the Bloch time that is 
required for the constant electric field to shift the kinetic 
momentum across the Brillouin zone (a ~ 0.142 nm is 
the carbon-carbon distance) and, of course, it is beyond 
the continuous Dirac model. The time scale At st appears 
for the tight-binding model as the time scale when the 
perturbation theory with respect to electric field breaks 
down (At s t ^> t 7 , where the microscopic time scale is 
< 7 = 7i/7 ~ 0.24fs, with 7 = 2.7 eV being the hopping 
energy) , and the dc response goes from a linear to a non- 
linear regime. 

The expression for (j 1 (t)) g given by (fTTf is a key for- 
mula in the study of the conductivity in the graphene at 
low carrier density beyond the linear response in dc. It 
describes the mean electric current of coherent carriers 
produced by the applied electric field. An exotic feature 
of the electronic transport at the charge neutrality point, 
as described by the Dirac model, is that one begins with- 
out any charges to be accelerated at all: there are no 
electrons or holes initially. The induced current can be 
considered, if one likes semiclassical style comments, as a 
consequence of two mechanisms: charged pairs of ultra- 
relativistic (with the constant velocity v F ) coherent elec- 
trons and holes are first created with a strong suppres- 
sion of large transverse momenta (\p 2 \ > \j e \E\ h/vp^j 
and after that their longitudinal kinetic momenta are co- 
herently increased by the electric field. The combined 
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effect of both processes that, in fact, cannot be localized 
and separated in the framework of QED, is described by 
Eqs. (|77p and ([78]) . In contrast to this, at high carrier 
density one has to consider the increment of the longi- 
tudinal kinetic momenta of an initially given number of 
incoherent carriers by the electric field. In the latter case 
the transport theory allows the semiclassical description 
of carriers. 

It follows from our results that the mean current in 
the graphene is parallel to the applied field, proportional 
to \E\ 3/2 , and grows linearly with time. The fact that 
the current grows indefinitely as T — > oo is a conse- 
quence of the absence of scattering and a backreaction 
mechanisms in the model of unlimited size under consid- 
eration: only effects caused by the applied external field 
have been taken into account. In the experimental situ- 
ation described in Ref. [19j, a constant voltage between 
two electrodes connected to the graphene was applied, 
and current-voltage characteristics (/ — V) are measured 
within ~ 1 s, which is a very large time-scale compar- 
ing with the ballistic flight time T ba i (the time which the 
electron spends to cross the material of a finite length 
L x ), Thai — L x /vp. To match our results with these con- 
ditions, our time T should be replaced by some typical 
time-scale that we call the effective time duration T e //- 
Some kind of dissipation process may truncate the pair 
creation at T^s, in which case T e ff = Tdis- The stan- 
dard candidates are collisions with impurities, phonons, 
ripplons, and the electron-electron interactions. A de- 
scription of how the pair creation is counterbalanced by 
dissipative processes so that a stable current settles down 
is still an open question. In the absence of the dissipa- 
tion, the transport is ballistic; in this case, considering 
a strip with lateral infinite width and a finite length L x , 
we assume the ballistic flight time T ba i to be the effec- 
tive time duration, T e // = T^ai- It is experimentally 
shown, see [H, [35[ and references therein, that the low 
temperature (9 ~ 5K) and gate voltage (\V g \ < 5 V) 
transport in current-annealed, suspended devices with a 
mobility ~ 170000 cm 2 /Vs is close to the ballistic limit 
over micron dimensions. In a realistic sample, placed on a 
substrate, the effective time duration T e ft can be many 
times smaller than T^ai, because of charged impurities 
or structural disorder of the substrate. However, such 
an effective time T e f j remains macroscopically large, so 
that Eq. ([75)) still holds. One can expect that in realis- 
tic high-quality (high mobility) samples the time T e ff is 
comparable with the time T ba i- In any case, for a finite 
flake length, the potential difference V = EL. X is finite, 
and one can consider the I— V of graphene devices. Then 
Eq. (|77p describes a regime where the current behaves as 

j ~ v 3 / 2 . 

This fractional power dependence of the current- 
voltage characteristic in graphene is at present called 
nonlinear (or superlinear) transport in graphene (I— V of 
the form j ~ V 3/2 ). In the Ref. 0, EE 111, this behav- 
ior was related to a possible pair creation by a constant 
electric field and experiments aiming at the observation 



of the effect were proposed. An experimental observation 
was recently reported It is agreed that such a su- 

perlinear transport is a distinctive feature characterizing 
the regime dominated by pair creation. 

It should be noted that our description of the quantum 
transport in graphene in the framework of strong-field 
QED is not restricted by a semiclassical approximation 
of carriers and it does not use any statistical assump- 
tions inherent in the Boltzmann transport theory. The 
estimations that follow show that this is important for 
the study of the conductivity close to the Dirac point. A 
typical density n^ r of carriers created in the ballistic case, 
given by Eq. ([76|) . e.g. for typical V ~ 1 V and L x ~ 1 
//m, is of the order n c g r ~ 6 x 10 11 cm~ 2 . Thus, in the 
general case, the density n c g r is of the same order as the 
impurity density n imp (10 10 -10 12 cm~ 2 [l3|). It should be 
stressed that this estimate is made for the ballistic case 
when T e ff = T bai . If T e// < T baU then the density n c g 
is less than rii mp . Moreover, the electric field creates the 
carriers in pure states with distribution (|14l) that differs 
significantly from the equilibrium distribution. The car- 
riers created at relatively small times T e ff ~ At st are 
wave packets with partial plane waves, which have wide 
range of wave lengths from oo to approximately 0.2 /jm. 
These lengths correspond to maximum kinetic momenta 
of the order \j e \E\ K/vp, see (fT4l . As T e f / increases the 
range of the longitudinal kinetic momenta of created par- 
ticle grows to its maximum e\E\ T e ff 3> \J e \E\ H/vf; at 
the same time the range of their transverse kinetic mo- 
menta remains unchanged. Thus, in the ballistic case, 
the lower bound of the range of the longitudinal wave- 
lengths decreases achieving its minimum of the order 4 
nm. The estimated minimum wavelength of the order 
0.2 fim is not much less than sample sizes that were used 
in experiments, and in some cases they are of the same 
order. Thus, we see that the Boltzmann theory, based on 
WKB approximation, is not well adapted to describe evo- 
lution of massless carries created from the vacuum. This 
is a nonstandard situation for usual transport problems 
in condensed matter physics. Our approach gives a non- 
perturbative description of the system evolution when 
quantum interference effects arc important. We see that 
the QED derivation of the superlinear I — V presented 
above gives additional arguments in favor of the interpre- 
tation that such a behavior is due to the pair creation. 
At the same time, in order to continue the study of the 
conductivity of graphene in the framework of the QED 
methods, one has to analyze in detail physical conditions 
in the graphene under which all the machinery, and es- 
pecially Eqs. ([77]) and ([75]). are valid. 



B. Mean electromagnetic field 

The mean time-dependent current (j 1 {t)) g is the 
source in the Maxwell equations for a mean electromag- 
netic field (E, B), where E = E + E ra£ ; and B = T5 ra d are 
electric and magnetic components, respectively. Here, 
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the initial external constant uniform electric field E sat- 
isfies the homogeneous Maxwell equations, and the fields 
E ra( i and ~B r ad are due to the current (j 1 (t)) g , repre- 
senting the backreaction of created pairs to the external 
field. The charged fermions in the graphene should feel 
the total field (E,B); in particular, the pair creation is 
induced by the mean electric field E. Since the current 
increases linearly as At increases, the backreac- 
tion to the external field should become relevant at some 
time instant. The situation looks similar to that studied 
in (loL [33j for QED in 3 + 1 dimensions, although there 
are considerable peculiarities in the case under consid- 
eration. Dealing with graphene devices, it is natural to 
assume that the constant strength E on the graphene 
plane is due to the applied fixed voltage V and there- 
fore is not changed when the created charges flow into 
the reservoirs which are located outside the graphene and 
have sufficiently large capacitances. It is assumed that an 
external current flows to the electrodes to maintain the 
fixed voltage, i.e., we are dealing with an open system of 
fermions interacting with classical electromagnetic field. 
Then Maxwell equations can be solved through infinite 
space and we may assume that the strength E is fixed 
by boundary condition at the infinity. The electromag- 
netic field is not confined to the graphene surface, z = 0, 
but rather propagates in the ambient 3 + 1 dimensional 
space-time, where z is the coordinate of axis normal to 
the graphene plane. To take this into account, we have to 
present the current term in the Maxwell equations as the 
current J ={J, 0, 0) restricted to a plane immersed in the 
3-dimensional space. We are interested in the time inter- 
vals At that satisfy the inequality < At < T. Then, 
taking into account that the current (j 1 (t)) g appears for 
At > 0, we can write J = (j 1 (*)>s 6(z)0(At) for At < T. 
The nontrivial Maxwell equations to be solved are: 

1 (9E $B 
VxS = , J + --, VxE = --, (81) 

where /j,q is the magnetic permeability and c is the speed 
of light in the vacuum. 

The only nonzero components of the irradiated field 
are the x-component of the electric field, E x ad , and the 
y-component of the magnetic field, By . Solving the 
Maxwell equations under initial condition E x ad = and 
B r y ad = at At = 0, one finds that for < At < T, the 
electromagnetic field produced by the current confined to 
the graphene sheet is: 

E rad = _^ sgn{E)A{cAt _ \ z \) 9 { C At - \Z\) , 

B™ d = S gn(z)E™ d /c. (82) 

The direction of the induced magnetic field changes 
sign across the graphene sheet, and sgn(z) = ±1 cor- 
responds to upper and lower regions. The electric field 
E x ad is opposite to the applied field E and is continuous 
in z-direction, in particular, at z = 0. A real graphene 
flake is a film of very large dimension in the x — y plane 



and a finite monolayer atomic thickness of approximately 
O.lnm to 0.2nm in the z-direction [13]. For realistic den- 
sities of carriers the potential field in z-direction looks 
like a deep potential well, which effectively forbids the 
motion of carries in this direction. Then it is natural to 
assume that the induced electromagnetic field inside the 
graphene is the mean between limiting values of the field 
from upper and lower regions as \z\ —¥ 0. Thus, the in- 
tensity of the induced electric field on the graphene plane 
reads: 

E r x ad = ~sgn (E) — \E\ 3/2 y^Jh At, (83) 

where a — /loce 2 /2h is the fine structure constant. As 
to the induced magnetic field, we believe that it is zero 
inside the graphene. 

The QED with an external constant electric field is a 
consistent model as long as the field produced by the in- 
duced current is negligible compared to the applied field, 
|-E , ™ d | <C \E\. This gives the consistency restriction: 

y/v F e\E\/HAt<: (84) 
4a 

In this case, the external electric field can be considered 
as a good approximation of the effective mean field. We 
call the typical time scale related to Eq. (|84|) . Atf,,. = 
Atstn/Aa, the time of backreaction. On the other hand, 
the dimensionless parameter in the l.h.s. of Eq. (|84|) sat- 
isfies the stabilization condition given by Eq. (f79j) . Thus, 
there is a window in the parameter range where the model 
is consistent, t 7 <§; At st <C Atfc r . Moreover, this restric- 
tion corresponds to a specific regime which might be rel- 
evant to some known experiments with graphene, as we 
shall see below. 

In the description of the carrier creation inside a 
graphene flake of a finite length L x , we use a simple pic- 
ture of the pair creation due to the T-constant field. We 
assume that the maximal duration of the electric field in 
our model is the effective time duration T e //, and that, 
in the case of ballistic transport, T e // = Thai- In typical 
experiments, L x ~ 1/im, and then T{, a i ~ 10 _12 s. Taking 
At = T ba i in Eqs. (JTHJ) and flBU), we obtain the following 
restrictions on the electric field: 

7 x 10 2 V/m< \E\ < 8 x 10 6 V/m. (85) 

Since the voltage is V = EL X , one finds the inequalities 

7 x 1(T 4 V < V < 8V. (86) 

These voltages are in the range typically used in experi- 
ments with graphene. 

In the experiment described in [19 ]. for instance, the 
I — V curve was studied in a range from decivolts to a few 
volts, in samples with lengths L x varying from 0.9 to 5.9 
fj,m and widths from 70 to 1500 nm. A power law / ~ V s 
was used to fit the I — V curves near the Dirac point, and 
it was found that l<<5<3/2. It was shown that the 
/ — V's in the graphene devices become superlinear in 
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the presence of disorder (in low mobility samples) while 
in high-quality (high mobility) samples, the superlinear- 
ity is masked by some other effect. The exponent S looks 
like a monotonically decreasing function of mobility for 
different devices and superlinearity vanishes for devices 
with high mobility, 5 — > 1. These results were interpreted 
as an interplay between pair creation and scattering by 
charged impurities and optical phonons. According to 
this interpretation, the superlinearity of the / — V in 
low mobility samples is in a qualitative agreement with 
expression (|77j) . while it is compensated by the contri- 
bution of the intraband current (the current of carriers 
that were present before the electric field was switched 
on, due to an imperfect experimental realization of the 
Dirac point conditions), which tends to saturate in high 
mobility samples due to interaction with optical phonons. 
The absence of a compensation of the superlinearity in 
low mobility samples is explained by the fact that the 
presence of a large number of charged impurities pre- 
vents the growth of the intraband current, making the 
interaction with optical phonons irrelevant. Except for 
this qualitative analysis, however, there is no theoretical 
description of these observations; the numerical analysis 
(in the NEGF approach) of 19J, in particular, does not 
include the optical phonons which are central to their 
interpretation of the effect, and it is argued that other 
effects should also be relevant for the analysis, as self- 
heating of the graphene sample, for instance |19| . 



In fact, according to our analysis, backreaction should 
not be neglected. We have shown that the consistency of 
the Dirac model with a given external field is ensured un- 
der the conditions ([79)) and (|84|) for At and E, which are 
thus extremely important for obtaining adequate physi- 
cal results in course of the calculations. Analyzing the 
experimental settings one can be sure that the stabiliza- 
tion condition (|79p is satisfied for all known measure- 
ments, which means that the lower bound in (|86p is re- 
spected. However, the condition (|84p could be violated 
because the voltage varies from 1 to 5 V, and sample 
length varies from 0.9 to 5.9 /im, as reported in [19j . 
Therefore, backreaction (appearance of the induced elec- 
tric field E^. ad ) cannot be ignored in calculating the mean 
current in high mobility samples. This means that, in the 
ballistic case, Eq. ([77)1 can give overestimated values for 
the current density of created particles for a given volt- 
age. In the case of low mobility samples, the effective 
duration time T e f f is due to the scale of dissipation pro- 
cesses Tdis, which can be many times less than the time 
Tbai, so that condition (|84|) holds true at At = Tdi S , and 
the regime of backreaction is not reached. As mobility is 
increased, backreaction becomes more pronounced, and it 
might compensate the superlinearity of the pair-creation 
process. In order to investigate whether this effect can 
lead to a transition from a superlinear to a linear I — V 
curve, we investigate in the next section how backreac- 
tion affects pair-production in graphene. 



C. Effective mean field and mean current 

In this section, we study a possible generalization of 
the above considered model with the T— constant exter- 
nal field. In particular, in the framework of this general- 
ization we are going to take into account the backreaction 
of the mean current to the applied electric field. 

Let us consider time intervals when the external elec- 
tric field is switched on, < At < T. In realistic cases, af- 
ter switching on, the external field E remains constant on 
the graphene plane due to the applied voltage V, which is 
supported by external sources to remain fixed. Accord- 
ingly, we assume that on the graphene plane E cannot 
be changed due to pair creation inside of the sample. As 
to the electric field E inside the graphene, we suppose 
that it can vary with time and it is directed along the 
axis x, so that E(i) = (E(t),0, 0). The x-component 
E (t) is a superposition of the applied external field E 
and a time-dependent electric field E^. ad (t) irradiated by 
the current induced in the sample, E(t) = E + E 1 x ad {t). 
In the initial time instant t\ — ti n , we have E(t\) = E. 
Let E{t) be a slowly varying field within a time interval 
Ati = ti + i — ti > such that the following condition 
holds: 



1 dE(t) 



E{t) dt 



A^«l, tefaU+i], (87) 



and let Ati be large enough to obey the stabilization 
condition of the type (IT9"|) . 



e \E(t)\ v F /hAti > 1 , te(ti,t i+ i] 



We consider now the time evolution of the mean field 
and the vacuum mean current between the initial time 
ti n when the external field is switched on and a time 
instant tfi n , which coincides with t ou t, or precedes it, 
tfin < t out . We first divide the interval (ti n ,tfi n ] into N 
equal intervals Ati, such that Ati = At2 = ••• = At at, 
J2iLi &U = tfin — tin- We suppose that Eqs. (j87J and 
(JHHJ) hold true for all intervals. That allows us to treat 
the electric field as approximately constant within each 
interval, E(t) E(ti), for t € 

Let us find the vacuum mean current density for the 
slowly varying electric field E (t) . We begin with the case 
t £ (iijia]- Here, according to Eqs. (|67j) and (|68|) . the 
current density (j 1 (t)) g is formed by the contributions 
from the vacuum states excited by the field E, these con- 
tributions having momenta that are restricted to the re- 
gion (|64[) . The phase volume of this region depends both 
on the magnitude of the electric field E and on the inter- 
val t— t\. It follows from the Eq. (|67p that the density 
(j 1 (t)) g grows as t — t% because the field E excites addi- 
tional states with larger longitudinal momentum p\. In 
the end of this time interval, when t = t2, the current 
density (j 1 (t)) g takes the form: 



Ah =sgn(E)D\E(t 1 )\ 3/2 At, 



(89) 
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with D = 2ir- 2 v 1 J 2 h- 3 / 2 e 5 / 2 . According to Eq. J83]), the 
field E™ d (t) has a direction opposite to E (ti) — E, and 
from inequality (|57|) its magnitude is small in comparison 
with it, \E r x ad {t)\ < \E(ii)\- Therefore, Eq. flU) deter- 
mines the leading term in the mean value (j 1 (t)) g in the 
corresponding large r-limit and in the mean field approxi- 
mation, and any corrections to (J (t)) g due to derivative 
d t E(t) are small. One can see that |£ , (t 2 )| < |-E"(ii)|. 
Now let us proceed to the case t £ (t 2 , £3] ■ The effective 
mean field E(t) is already different from E and, as before, 
the field E(t) is approximately constant, E(t) w E(t 2 ). 

The constant field E{t 2 ) acting on the interval (£2,^3] 
excites new vacuum states which do not participate in 
the formation of the current density (j 1 (t)) g on the pre- 
vious step. These new states have longitudinal momenta 
PiSgn(£') greater than states excited from the vacuum 
during the time £2 — t\ and transversal momenta p 2 lim- 
ited by the field E(t 2 ). The region of these momenta is 
defined by inequalities similar to those in Such in- 

equalities in the mean field approximation have the form 

1/2 



\p 2 \ < Je\E{U)\h/v F 



e\E(U)\v F /h(t-U)- K 



U + 



K 



\eE{U)v F /h 



1/2 



< -- 



Pi 



eE{U 



< t - 



K 



\eE{U)v F /h 



1/2 ' 
(90) 



In the case t £ (^2,^3], one has to set ti = t 2 in these 
relations. As a result, the current density (j 1 (t)) g in the 
second interval takes the form 



0" 1 (*)> ff = Aji + Aj(t - 1 2 ) . 



Aj(t - t 2 ) = sgn(£)D \E(t 2 )\ 3/2 (t - t 2 ) 



(91) 



in the large r-limit. As in the case considered above, this 
current irradiates the corresponding field E™ d , which is 
directed against to the constant field E(t 2 ), but its mag- 
nitude is less than the latter. Because t — t 2 < At 2 and 
it is supposed that condition (1571) holds true, one can 
see that the contribution from the derivative d t E(t) to 
the current density (j 1 (i)) 9 is much less than the leading 
term, given by (liTTl) . Due to the contribution E x , the 
mean field will decrease, \E(t)\ < \E(t 2 )\, reaching its 



minimal value at t = t^. The further evolution of the 
vacuum mean current density and mean field has similar 
behavior for time t from any interval Ati. In each such 
intervals the region of the momenta is given by Eq. (|90p . 
In the general case (in the corresponding large r-limit) 
when t E (tMitM+i], 2 < M < N, we can represent the 
current density (j 1 (t)) g as the following sum of partial 
contributions: 



M-l 



(i 1 (*)>s = Aii + Ai(t-t M ), 



1 3/2 



At, 



Aji =sgn(E)D\E(U 

Aj(i - t M ) = sgn(E)D \E(t M )\ V2 (t - t M ) . (92) 

Condition (|87|) guarantees that corrections to the lead- 
ing terms Aji and Aj(t — tu), given by (|92|) . due to the 
derivatives dtE(t) are small in each interval (ti,ti + {\. In 
each time interval, the mean current irradiates the corre- 
sponding field E r x ad directed against the external field E, 
thus, the mean field E(t) is a monotonically decreasing 
function for t 6 (ti n ,tfi n ]. 

Let us use the approximations 



\E(U)\ V2 AU 



\E(tM)\ ^ (t — tM) 



\E(t)\ 3/2 dt, 



f \E(t)\ 3/2 dt. 



Then the mean value (j (t)) g can be approximated by 
the following integral form: 



(j\t)) g = S gn(E)D 



\E{i)f 2 di. 



(93) 



tin 



Note that Eq. (|93|) holds true within large r-limit, in 
which we neglect a contribution from the derivative 
dtE(t), if the condition (|57J) is valid for all the intervals 
(ti, U+i]- 

In a similar manner, one can find a unified represen- 
tation of the vacuum mean current density (j 1 ^))^ and 
the diagonal elements (T MM (i))g for the slowly varying 
electric field E (t) in the same approximation: 



(j\t)) g =i8eP 1 AV(x,x')\ x=x , , (T^(t)y g =i8P^(x,x')\ 
A p (x,x) — —i j h±(x±,x'_ L )h^(xpx'i [ \)dt , 



(94) 



t in +KAt s 



h \\( x b x '\\) 





E(t) 




\Mt) 




Mt)\ 



exp 



-eA 1 (t){x 1 -x' x ) 



h ± (x ± ,x x )=[ w ) (27r)2 exp 



(x 2 -x' 2 ) 2 e\E(t)\ 



(95) 
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where 

-m^ 1 =e[\Mt)\-\A 1 (t)\] , 
vf at li i i u 



Ai{t) = J E(t)di+EU 



1/2 

At t = (e |-E(t)| Vjr/h) , and Ai st is the characteristic 
time-interval determined by Eq. (|5U|) . These equations 
are a generalization of representation (l68l) for the case 
under consideration, where A p (x, x') is a generalization 
of the function A p (x, x') in (|67l) . One can see that the ex- 
pression (j 1 (t)) g given by (IMl) coincides with that given 
by Eq. (1931) in the large r- limit. In this approximation 
the mean values (T uu (t)) g have the form: 

(Too (t)) g = (Ti 1 (t)) g = (T u (t))> 1 
(T 22 {t)) g = (T 22 {t)Y g + (T 22 {t)Y g , 



with 



pt-K At t 

(T n (t)) p g =2ev F p \E(t) 

Jt in +KAt st 



3/2 



x [\At(t)\ - dt 

r t-KAi t 



(T 22 (t))P =^hp 



1 





E(t) 


5/2 


1^1 (t) 




^l(*)| 



(T 22 (t)) c = ~£ R (3/2) hp\E (t) 



andp = vr- 2 « F 1/2 ft- 3 / 2 e 3 / 2 . 



1 3/2 



(96) 



By using Eq. (1931) , we obtain that in the Maxwell equa- 
tions for the mean electromagnetic field the current term 
for t < tfin has the form J = (j 1 (t)) g S(z)6(t - h). The 
only nonzero components of the radiated field are the x- 
component E x ad of the electric field and the y-component 
By ad of the magnetic field. Solving the Maxwell equa- 
tions (|8T|) under the initial conditions E x ad — and 
B r y ad = at t = t in , one finds that for < t - t in < 
tfin— Urn an d within a distance \z\ smaller than c(t — t in ) 
to the plane, the electromagnetic field produced by the 
current confined to the graphene sheet is: 



ET(t,z) 



sgn(E) 
x B(t- 



p D 



E{t-l/c)\ 3/2 
l/c)6(l- \z\)dl, 



B: ad (t,z) =sgn(z)E r x ad (t,z)/c. 



(97) 



It can be seen that the radiated field is a function of only 
the light-cone variable, E r x ad (t, z) = E r x ad (t ~ \z\ /c). 

The electric field E x ad (t, z) has the same limit value 
near the plane from the upper and lower regions and we 
assume that E x ad (t,z) is continuous at z = 0. Thus, 
inside the graphene film the irradiated electric field is 
E r x ad (t) = E r x ad (t, 0). The mean magnetic field inside the 



graphene film is zero. Extracting the derivative dtE(t) = 
d t E r x ad (t) from flE), we get that for t G (t in ,tf in \ a 
self-consistency equation for the mean field inside the 
graphene has the form: 



d\E(t)\ 
Jt 



3/2 



t G (U 



(98) 



As was already established, the mean field field E(t) is 
collinear to the external field E. Taking into account the 
initial condition, E(ti n ) — E, one finds the solution of 
equation (|9"5)l as follows: 



E(t) = Ee(ty 



<t) = 



2a , . , . 

1 + — (t-t in )/At s 



(99) 

Then the irradiated field is E r x ad (t) = E(t) - E. _ 

As t — > tfi n the magnitude of the mean field E(t) de- 
creases. Therefore, condition (|88|) will be satisfied for all 
the intervals (tj,fcj+i] if it is satisfied for the last interval 
(t ff, tfin]- in our approximation, this condition defines 
the minimal possible length of the intervals Ati, 



Ati > At = At st e(t f i 



(100) 



Condition (|57|) is satisfied for all the intervals (ti,t i+ i], 
if it is satisfied for the period (ti,t 2 ]. With Ati 3> At, 
the latter implies the following restriction for the interval 



t fin 



4a 



(tfin) -C 1 => {tfin — U n )/At 3t <C 



(101) 



The mean field approximation is consistent if inequality 
(| 1 1 [) holds true. This restriction is much weaker than 
the one given by Eq. (|84p. which ensures consistency of 
the external field approximation. Thus, in the approach 
under consideration the largest time scale is At fi n , which 
is related to time scales considered before as follows: 



At 



fin 



7T 71 

Atbr = —^;At s t . 



2a 



When tfi n ~ tin approaches At^„, the magnitude 
|-E(f/i„) decreases and becomes of the order |-E m i n | = 
\E\ (4a/7r) 2 . One can verify that in the case of ballistic 
transport, for all the typical experimental parameters, 
when the voltage varies from 1 to 5 V and sample lengths 
vary from 0.9 to 5.9 pm (see [19]), the corresponding bal- 
listic times satisfy condition (jlOip . T^ai <C Atf in . This 
means that in realistic cases the evolution of quantum 
states in course of the external field action satisfies the 
restrictions which justify the approximation under con- 
sideration, so that we can set t fi„ = t out . 

Note that in the case when the field duration t — t in ex- 
ceeds At fi n (that is unattainable for samples available at 
present), the mean field E(t) becomes too weak in com- 
parison with the external field E, and then the model 
with pair creation by a constant external field fails to 
work. It is possible to consider the evolution of the mean 
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field and the vacuum mean current for time intervals 
greater than Atf in if one is able to calculate explicitly 
pair creation from vacuum due to a time-dependent elec- 
tric field, which could be a task for further study. How- 
ever, it is natural to assume that the residual effect of a 
weak electric field cannot significantly change the asymp- 
totic behavior achieved by the time t = tti n , which ef- 
fectively means that the effective external field duration 
is T ss tfi n —ti n . That is why in what follows we set 
tfin = tout for any t out . 

Substituting Eq. ((Ml) in Eq. O, we find the following 
result for the density of the vacuum mean current: 

O' 1 ^)}* = ^ [i - <t)~ 2 } » te(t in ,t fin }. (102) 

Thus, we have found self-consistent solutions for mean 
field and vacuum mean current. From Eqs. and (|102[) 
we see that at (t — ti) 3> Att, r = A£ st 7r/4a, the density 
of the vacuum mean current and radiated electric field 
take asymptotic forms: 



(i x (i)) £ 



2nha ' 



E™«(t)n-E. 



(103) 



Thus, the self-consistent system of the mean field and the 
vacuum mean current at a given external electric field E 
tends to a dynamic equilibrium state in which the exter- 
nal field inside the graphene is completely compensated 
by the radiated electric field. In this state the parti- 
cle production is stopped and the vacuum mean current 
saturates. Close to this regime, the I — V is almost lin- 
ear. The momentum transferred by the applied force is 
limited by the finite value e \E\ At st ~K jloi. Due to the 
discrete structure of the lattice, formed by atoms at the 
carbon-carbon distance a, there is a natural momentum 
cutoff of the order H/a. For this reason, instead of being 
restricted by the Bloch time At b mentioned in subsection 
IIV A[ the applicability of the Dirac model is restricted, 
because of the crystalline structure of graphene, by the 
field strength Eb, 



— — - 3 x 10 7 V/m. 
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Substituting Eq. (|99j) into Eq. (|96|) . we can explicitly 
find the diagonal elements (T uu {t)) p g . For example, the 
mean energy density reads: 



1 



1 



(Too(t)) g = 8ev FP \E\^(At br f [3 ^ , ^ 

(104) 

For (t — t in ) ^> Atb r , we obtain from (|104[) the following 
asymptotic form: 

(Too (t)) g K~ev F p\E\ 5/2 (At br ) 2 . 

One can see that (Too(i)) 9 ~ |i?| 3 ^ 2 asymptotically for 
large times, while for small intervals (t — tj„) <C Atbr one 



has a dependence of the form (T 00 (t)) g ~ |i?| 5 ^ 2 . The 
asymptotic behavior of the element (T22{t)) p g reads: 



hp 
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The behavior of the term (T22(t)) c g , given by (l96)l . is com- 
pletely determined by the field 2£(i) at a given time in- 
stant, which means that this term tends to zero asymp- 
totically as e(t)~ 3 , and therefore can be neglected. We 
see that (T 00 (t)) g = (T u (t)) fl > (T 22 {t)) g during the 
whole evolution up to asymptotically large values of the 
interval (t — ti n ). Then, in this approximation, we have 
that trace (Tjf(t)) g = 0, as it can be expected for the 
equation of state for massless particles. 

At t = tout, the expressions (|102|) and (I104j) represent 
the current density and energy density of created parti- 
cles, respectively, and the consistent mean electric field in 
graphene is given by Eq. (1991) . According to our problem 
setting, before the time moment t outl we have an unitary 
evolution of a pure state of the Dirac-Maxwell system. At 
the initial time instant tin, this pure state is the vacuum 
for carriers and the coherent state of electromagnetic field 
with the initial mean value E. In the strong-field QED, 
"the measurement" which produces the decoherence (as 
a result of which we obtain many-particle state of car- 
riers and a final state of the electromagnetic field speci- 
fied by the mean value E (t out )) occurs not necessary at 
the moment when the external field switches off. This 
measurement can be done in any time instant after t out 
because the further evolution of the system is trivial and 
mean values remain unchanged. In the case of the finite 
graphene size and fixed constant voltage this is not true 
because in such a case we adapt our model with the T- 
constant external field to the situation where the effective 
duration is given, i.e., we set T = T e ff and, therefore, by 
definition, we have t out = U n + T. (Note that one has 
to know the times i,- n and T in a different experimen- 
tal situation for the measurement of the time evolution 
of the vacuum mean current, which has been been early 
proposed [4| but not implemented yet.) In this case, the 
time instant t out of switching off the effective electric field 
is understood as the effective time instant of the decoher- 
ence. We assume that the decoherence stops evolution of 
pure states in a sufficiently short time, so that the final 
mean electric current and EMT do not feel the effect of 
switching off the effective electric field. 

When the effective duration T e ff is due to the time 
scale of dissipation process Tdi S = l m fp/ v F, we set t out — 
tin = T e ff = lmfp/vF, while in the ballistic case l m f p — 
L x and T e ff — T^ai- Then, using (|102|) . we obtain the 
I — V curve in the following form: 




(105) 
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Example 1: Considering a sample with L x ~ 1/xm and 
a voltage IV in the ballistic case, we have T^i — 10~ 12 s 
and \E\ — 10 6 V/m. Alternatively, in the case of dissipa- 
tion, we can consider a sample of arbitrary length, where 
the effective duration is T e ff — 10~ 12 s, and the applied 
field strength is \E\ = \V\ /L x = 10 6 V/m. In both cases, 
2aT e ff(irAt st )~ 1 = 0.18, and condition flM]) is satisfied. 
But then backreaction can be neglected and the I — V 
curve from eq. (| 1051) is close to j ~ V 3 / 2 . 

Example 2: Considering a sample with L x ~ 4/im 
and a voltage 4V in the ballistic case, we have 
2aT eff (nAt st )- 1 = 0.72. Thus, condition flM]) is not 
satisfied, and then the backreaction contribution in eq. 
(|105l) is important, resulting in an almost linear I — V. 
We see that our approach describe a transition from a 
superlinear to a linear I — V curve when passing from 
low-mobility to high-mobility samples, in a manner sim- 
ilar to that observed in 19]. 

Note that outside the graphene film there is an irra- 
diated electromagnetic plane wave E x ad (t, z) = E x ad (t — 
\z\/c) of linear polarization given explicitly by Eqs. (|97|) 
and (|9"9")l . Therefore, the energy dissipation due to this 
irradiation must be taken into account. In principle, such 
irradiation may be experimentally observed. It is the ra- 
diation due to the time-dependence of the mean current. 
We have calculated it nonperturbatively in the case when 
its electric component strength is comparable to the ex- 
ternal field strength. Thus, the mechanism of the energy 
dissipation due to irradiation considered here differs es- 
sentially from that discussed in Ref. 18]. The authors 
of the latter work have treated the radiation due to the 
electron-hole recombination in the lowest order of the 
perturbation theory assuming large enough densities of 
carriers. That is why the relevance of such a considera- 
tion to the quantum transport close to the Dirac point is 
not clear. 



V. SUMMARY 

In this summary, we briefly list the main new results 
obtained in the article and add some relevant comments. 
These results are collected in three blocks: I-General re- 
sults in strong-field QED; II- Adaptation of the general 
results to the Dirac model in the graphene; III- Analysis 
of some immediate consequences to the graphene physics. 

I. General results in strong-field QED 

The one-loop renormalized mean current and EMT in- 
duced in the vacuum by the T-constant external electric 
field are computed in the framework of strong-field QED 
in spaces of arbitrary dimensions. 

In the large r- limit, these quantities are represented 
as sums of local contributions due to the vacuum po- 
larization, {{j p {t)) c L and Re (r MI ,(£))£ e „), and of global 
contributions due to the vacuum instability, (Re 



and Re (T^ (*)}*>), 

(f{t))=(j»{t))±+Re(j»{t))\, 
<V(t)) re „ = Re (T^(t)) c ren + Re {T^{t)) p ■ 

These contributions are studied in detail. The vacuum 
polarization contributions to the EMT arc expressed via 
the real part of the one-loop effective Euler-Heisenberg 
Lagrangian. In odd dimensions, unusual peculiarities of 
the vacuum polarization emerge: along with the lon- 
gitudinal mean current, Re(j^(i))j| (which behaves in 
the same way as in even dimensions), there appears a 
transversal mean current, (J^{t))±. Its components cor- 
respond to fermions of different chiralities (pseudo-spin) 
moving in opposite directions. The sign of this current 
depends on the fermion species. That is an indication 
that the Chern-Simons term is present in a properly 
regularized effective action with odd number of fermion 
species. In the case of an external electric field, this trans- 
verse mean current depends essentially on the fermion 
mass (in contrast to the corresponding contributions to 
probability amplitudes of processes); in particular, it is 
zero in the massless case. Thus, the vacuum polarization 
in the electric field in odd dimensions is qualitatively dif- 
ferent for mean values and for amplitudes of processes. 
This is the reason why the case of electric-like field differs 
substantially from the case of magnetic-like field. This 
fact is important in understanding properties of new ma- 
terials where the electronic structure is described by the 
Dirac model of (2+l)-d fermions, in particular in the case 
of the interface transport in topological insulators, where 
there is only one fermion species. Moreover, it should be 
noted that even in the case of even number of fermion 
species with nonzero masses, neutral electron-hole pairs 
induced due to the applied electric field will be concen- 
trated near lateral surfaces. Pairs with different chiral- 
ity are concentrated on opposite lateral surfaces (there 
is a polarization with respect to the chirality). Other- 
wise, properties of these mean values in odd and even di- 
mensions do not differ essentially. In particular, the lon- 
gitudinal mean current components Re(j^(i))jj are lin- 
ear functions of T, and the energy density (T 00 (t)} ren 
and the pressure along the direction of the electric field 
(Tn(t)} ren are quadratic functions of T. 

II- Adaptation of the general results to the Dirac model 
in the graphene. 

The general results described above are used in the 
non-perturbative consideration of electronic and energy 
quantum transport in graphene at low carrier density and 
low temperatures. 

We have found: the time-dependence of the mean 
longitudinal current density (j 1 (t)) g , the mean EMT 
{Tfifj, (t)) g , and the total number density of electron-hole 
pairs n C g created by the electric field in the graphene. 
In these calculations, we have used the strong-field ap- 
proximation, the large r-limit (when the duration T of 
the external electric field is sufficiently large), and the 
backreaction of created particles in screening the exter- 
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nal field was taken into account. All these quantities are 
obtained for the case of an unitary evolution of a pure 
state of the Dirac-Maxwell system. To adopt our consid- 
eration to the existing experimental situations, where the 
length of graphcnc flakes is finite and a constant voltage 
is applied, the time T was replaced by some effective time 
T e ff, which is determined by a certain decoherence. As- 
suming that this decoherence stops the evolution of pure 
states in sufficiently short times one can neglect effects of 
switching off the effective electric field on the final mean 
electric current and EMT. 

The new approach proposed for treating the quantum 
transport in graphene allows to study unitary (coherent) 
evolution of system in the regime of large duration of a 
strong external field, when the electrodynamic backre- 
action is essential. In spite of the fact that we do not 
consider statistical effects of the decoherence, we believe 
that exact analysis of results of unitary evolution will 
be helpful as an initial stage for further considering the 
thermalization of the system. 

It should be noted that in some high-energy physics 
problems (QED and QCD) in (1 + 1)- and (3 + 1)- 
dimensions the backreaction from the pair production 
in electric fields was frequently discussed, in particular, 
in the framework of WKB approximation (see e.g. 1361 - 
and using the numerical approach (see e.g. [39j|). 
Recently, the exact solution for massless QED in 1 + 1- 
dimensions was obtained, which by the definition takes 
the back reaction into account completely [4(|. The au- 
thors of all these works considered only closed systems of 
fermions interacting with a mean electromagnetic field, 
in which the total energy was conserved. However, in 
considering the graphene in the framework of the Dirac 
model the backreaction due to the pair creation from the 
vacuum has never been studied. Our consideration of 
the backreaction in the graphene due to the pair creation 
is completely new and, moreover, is essentially different 
from that of the works cited above, because it concerns 
an open system with a source that maintains a given 
voltage. 

Ill-Analysis of some immediate consequences to the 
graphene physics. 

In the large r- limit, the mean values (j 1 {t)) g and 
(T^j {t)) g are proportional to the number density of pairs 
r c J of positive and negative charged states excited by 
the electric field per unit time. At the end of pure 
state evolution these quantities represent the current den- 
sity and EMT of created particles, respectively, so that 
{j\t)) g ~ \Ef 2 and (Too {t)) g = <T U [t)) g ~ |£| 5/2 . 

It is shown that there exists a parameter range 1 <C 
T/At st <C n/Aa for which the external electric field of 



long duration T is a good approximation of the effec- 
tive mean field inside the graphene. Two corresponding 
characteristic time scales are established: one of them 
is At st = (vfs\E\ /h) , the time duration when the 
dc response goes from a linear to a non-linear regime 
and the effect of the real electric field is indistinguishable 
from the effect produced by a constant field; the other 
one is Atb r = Ai s t7r/4a, the time duration when the 
backreaction becomes important. Such a specific regime 
is relevant to describe the experimentally observed super- 
linearity j ~ V 3 / 2 of the I — V in low mobility samples. 

A generalization of the Dirac model with the In- 
constant electric field is constructed that takes into ac- 
count the backreaction of the mean current to the ap- 
plied electric field set by a constant voltage. It is shown 
that in the case of graphene the electrodynamic backre- 
action is relevant and one can see the interplay between 
two and three dimensions: the time-dependent mean cur- 
rent is confined to the plane, but its radiation, form- 
ing the backreaction on the plane, escapes to the three- 
dimensional space in the form of linearly polarized plane 
electromagnetic waves. A self-consistent solution of the 
Dirac-Maxwell set of equations for this generalized model 
is found and the effective mean field and effective mean 
values of the current and energy-momentum tensor are 
calculated. In this case the self-consistent system of the 
mean field and the vacuum mean current tends to a dy- 
namic equilibrium state in which the external field inside 
the graphene is completely compensated by the radiated 
electric field, the particle production is stopped, and the 
vacuum mean current saturates. Close to this regime, the 
I — V is almost linear and the mean energy density reads: 
(Too (t))g ~ |i?| 3 ^ 2 . This mean field approximation is 
consistent if the inequality T/At st <C n 2 /8a 2 holds true. 
This restriction is much weaker than the one that deter- 
mines consistency of the external-field approximation. 

We show that non-linear and linear I — V experimen- 
tally observed in low and high mobility samples, respec- 
tively, can be explained in the framework of the presented 
consideration and that such a behavior is a consequence 
of the fact that the conductivity in the graphene is es- 
sentially due to the pair creation from vacuum by the 
applied electric field. 
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